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^ ' Abstract 

> ; 

O , We estimate the accuracy with which the coefficient of the CP even dimension six operators involving 

. Higgs and two vector bosons {HVV) can be measured at hnear e"'"e~ coUiders. Using the optimal 

. observables method for the kinematic distributions, our analysis is based on the five different processes. 

First is the WW fusion process in the t-channel (e"'"e~ — > Ve^eH)^ where we use the rapidity y and the 

' transverse momentum px of the Higgs boson as observables. Second is the ZH pair production process 

Cii' 

I ' in the s-channel, where we use the scattering angle of the Z and the Z decay angular distributions, 

Oh' 

I reproducing the results of the previous studies. Third is the t-channel ZZ, fusion processes (e^e~ — > 

e'^e~H), where we use the energy and angular distributions of the tagged e"*" and e~. In the fourth, we 
[ consider the rapidity distribution of the untagged e^e~H events, which can be approximated well as 

. the 77 fusion of the bremsstrahlung photons from e+ and e~ beams. As the last process, we consider 

. the single tagged e'^e~ H events, which probe the 76^ He^ process. All the results are presented 

■ in such a way that statistical errors of the constraints on the effective couplings and their correlations 

00 ' are read off when all of them are allowed to vary simultaneously, for each of the above processes, 

^ ; for = 120 GeV, at ^ = 250 GeV, 350 GeV 500 GeV and 1 TeV, with and without e" beam 

^P, I polarization of 80%. We find for instance that the HZZ and HWW couplings can be measured with 

\ 0.6% and 0.9% accuracy, respectively, for the integrated luminosity of L = 100 fb~^ at ^/s = 250 GeV, 

^ ! 350 GeV and L = 500 fb~^ at ^/s = 500 GeV, 1 TeV, for the luminosity uncertainty of 1% at each 

energy. We find that the luminosity uncertainty affects only one combination of the non-standard 
couplings which are proportional to the standard HWW and HZZ couplings, while it does not affect 
the errors of the other independent combinations of the couplings. As a consequence, we observe 
that a few combinations of the eight dimension six operators can be constrained as accurately as the 
two operators which have been constrained by the precision measurements of the Z and W boson 
properties. 
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I. INTRODUCTION 



The Standard Model (SM) of the elementary particles based on the SU(3) x SU(2) x 
U(l) gauge symmetry has proved to be a successful theory to interpret all the precision data 
available to date. SM predicts a light Higgs boson whose discovery is one of the prime tasks of 
the upcoming future colliders. 

In fact, the present electroweak precision measurements indicate the existence of a light 
Higgs boson [l|, S]. Experiments at the CERN Large Electron Positron collider (LEP) set the 
lower bound on its mass of 114.1 GeV at the 95% confidence level (CL) [l|. The Fermilab 
Tevatron, which collides proton and anti-proton at y/s = 2 TeV, is currently the only collider 
which can produce low mass Higgs bosons. Analysis with Run lib data samples by the CDF 
and D0 detectors indicates that the Tevatron experiments can observe the Higgs boson with 
about 10 fb~^ total integrated luminosity for the mass of around 120 GeV 3|. The Large 
Hadron Collider (LHC) at CERN will start colliding two protons at y/s = 14 TeV in the year 
2008, and is geared to detect the Higgs boson in gluon-gluon and vector-boson fusion processes. 
It will measure ratios of various Higgs boson couplings through variety of decay channels at 
accuracies of order 10 to 15% with 100 fb^^ luminosity 4|. 

Despite the success, SM presents the naturalness problem due to the quadratic sensitivity 
of the Higgs boson mass to the new physics scale at high energies, which implies that there is 
a need of subtle fine tuning to keep the electroweak symmetry breaking theory below the TeV 
scale. To put it in another way, this may suggest an existence of a new physics scale A not far 
above the TeV scale. The key to probe the new physics beyond the SM theory is to clarify the 
origin of the electroweak symmetry breaking, the Higgs mechanism. Therefore, it is necessary 
to measure the Higgs boson properties as precisely as possible, especially the HVV couplings, 
because they are expected to be sensitive to the symmetry breaking physics that gives rise to 
the weak boson masses. 

With this motivation, we re-examine the potential of the future e+e~ linear collider, the 
International Linear Collider (ILC) in the precise measurement of the HVV couplings. Clean 
experimental environment, well defined initial state, tunable energy, and beam polarization 
renders ILC to be the best machine to study the Higgs boson properties with high precision. In 
this paper, we study the sensitivity of the ILC measurements on all the HVV {HWW, HZZ, 
HZ'-f and H'-f-f) couplings comprehensively and semi- quantitatively by using all the available 
processes with a light Higgs boson {tuh^ 120 GeV); e~ Ue^eH with t-channel W ex- 
change, e"*" e" ^ Z H with s-channel Z exchange, e"*" e~ e~ H with t-channel Z exchange. 
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no-tag e"*" e~ — (e"*" e~) if process from 77 fusion, and single-tagged e"*" e~ — > (e^) if process 
that probes ^ H via t-channel 7 and Z exchange. 

In order to quantify the ILC sensitivity to measure various HVV couphiigs simultaneously, 
we adopt the powerful technique of the optimal observables method It allows us 

to measure several couplings simultaneously as long as the non-standard couplings give rise to 
different observable kinematic distributions. The results can be summarized in terms of the 
covariance matrix of the measurement errors, from each process at each energy, that scales 
inversely as the integrated luminosity. 

In order to combine results from different processes and at different energies, we adopt the 
effective Lagrangian of the SM particles with operators of mass dimension six to parametrize 



all the HVV couplings 



This allows us not only to compare the significance of the 



measurements of various HVV couplings at different energies and at different colliders, but also 
to study what ILC can add to the precision measurements of the Z and W boson properties in 
the search for new physics via quantum effects. We therefore parametrize the HVV couplings 
as linear combination of all the dimension six operators that are allowed by the electroweak 
gauge symmetry and CP invariance. 

Some of the previous studies based on the optimal observables method are found for CP- 



violating effects m e~ ^ Z H via HZZ and HZ'-j couplings |8|, ll2|, and also in e"*" e~ ttH 
8|. CP conserving and CP violating effects in e"*" e~ — ZH process has been studied in 
ref. jl3, [1^. In refs. [s], Is], 1^ all the relevant couplings are varied simultaneously, and their 
correlations are studied. More recently, the ILC sensitivity to the HZZ and HWW couplings 
has been studied in refs. Is], [igI, Bounds on the coefficients of the Higgs- vector boson 
dimension-6 operators have been found in refs. [isl, [l^ based on non-observation of the Higgs 
boson signal at the Tevatron. Whenever relevant, we compare our results with the previous 
observations. 

This paper is organized as follows. In section [Til "we describe the low energy effective in- 
teractions among the Higgs boson and the electroweak gauge bosons arising from new physics 
that is parametrized in terms of the effective Lagrangian of the SM particles with operators 
up to mass dimension six. In section IIIIl we introduce the optimal observables method and 
explain how we perform the phase space integration when some of the kinematic distributions 
are unobservable, such as neutrino momenta and a distinction between quark and anti-quark 
jets. Although we present numerical results for unpolarized beams and for 80% polarized e~ 
beam only, all the formulas are presented for an arbitrary polarization of e~ and beams. 
After introducing final state cuts, such as those for the tagging and those for selecting or ex- 
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eluding Z ^ f f events, we present the total cross sections for all the five processes at ^/s=200 
GeV-1 TeV for =120 GeV, and at ^=250 GeV, 500 GeV, 1 TeV for = 100-200 GeV. 
Then in section HV] we compute the statistical errors of the non-standard HWW couplings 
extracted from measurements of the ly VT-fusion process, e^e^ i^eH. In section |Vl we 

study the constraints on the HZZ and HZ'-y couplings extracted from ZH production. In 
section IVll not only the HZZ and HZ^ couplings but also the if77 coupling are studied in the 
double-tag e^e~ — e^e~ H process via t-channel Z and 7 exchange. In section IVIIt we obtain 
the constraints on the ^^77 coupling from the 77 fusion, in no-tag e^e~ — > {e~^e~)H events, 
using the equivalent real photon approximation. In section IVIIIt we consider the single-tag 
e'^e~ (e^) H process to constrain the H'-ff and H^Z couplings. In section [IX| we address 
the implication of luminosity uncertainty on the measurement of these couplings. In section |Xl 
we summarize all our results, compare them with previous studies, and present our estimates 
for the ILC constraints on the dimension six operators, which are then compared with the 
constraints from the precision electroweak measurements of the W and Z boson properties. In 
Appendices we present our parameterizations of the 3-body phase space (Appendix |A]) , and 
the explicit forms the t-channel and s-channel currents and their contractions that appear in 
the helicity amplitudes (Appendix [B]) . 



II. GENERALIZED HVV VERTEX WITH DIMENSION SIX OPERATORS 

In our study, we adopt the effective Lagrangian of the Higgs and the gauge bosons with 
operators up to mass dimension six, 

i^eff = LSM + 5^ \ (1) 

i 

Where Lsm denotes the renormalizable SM Lagrangian and ^'s are the gauge-invariant op- 
erators of mass dimension 6. The index i runs over all operators of the given mass dimension. 
The mass scale is set by A, and the coefficients f^^^ are dimensionless parameters, which are 
determined once the full theory is known. Excluding the dimension 5 operators for the neutrino 
Majorana masses, and the dimension 6 operators with quark and lepton fields, we are left with 
the following eight CP even operators that affect the HVV couplings. Notation of the operators 
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are taken from the reference [2 



Oww 
Obb 

Ow 
Ob 
04,1 



(<i>t$)(D^<i.)t(D'^$), 

^(9^(<l>'^<l>)9^(<l>"^$). 



(2a) 
(2b) 
(2c) 
(2d) 
(2e) 
(2f) 

(2g) 
(2h) 



|, and the covariant derivative 



Here $ denotes the Higgs doublet field with the hyper-charge Y 
is = + igwT°'W^ + igYYB^, where the gauge couplings and the gauge fields with a 
caret represent those of the SM, in the absence of higher dimensional operators. The gauge- 
covariant and invariant tensors W ^i, and B^y, respectively, are W ^i, = igwT"'W^^,, and B^^, = 
^gyBfj^iy. The coefficients of the operators ( l2al) -( l2hl) . which are denoted as ff'^/h? in the effective 
Lagrangian of eq.([T]), should give us information about physics beyond the SM. So far, the 
precision measurements of the weak boson properties constrained the operaotrs Obw and 
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Oaii, which have been useful in testing some models of the electroweak symmetry breakdown 
23|. In this report, we explore the accuracy with which the ILC experiments can measure the 
coefficients of all these eight operators when a light Higgs boson exists. 

When the Higgs field acquires the vacuum expectation value < $ >= :^(0, f )^, the bilinear 
part of the effective Lagrangian of eq.([T]) is expressed as 



^eff 



; 1 + 



2 A2 



WW 



yra jyafiu 



2 A2 



BB 



+^ 1 + 



2 A2 



(Li + U {gw W\ - gy B,) {gw - gv B^ 



+ ^ 1 + 



(/0i + /04 + 2/^2)) [d,Ho) (to 



2 A2 

After renormalization of gauge fields and their couplings. 



2 H? 







gw_f_ 

4 A2 



/l 



WW 



gl v"^ 

4 A2 



BB 



B. 



At , 



gw w; = gw w;-, 

gv B^ = gy Bfj_; 



(3) 

(4a) 
(4b) 
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and after diagonalization of the mass squared matrices, the effective Lagrangian reads 



^eff 



_ -w^ W 

2 /j'f 



1 y ytiv 



2 ^ 2 



ml + 



(5) 



where 



= 


9w V' 
4 


m| = 


9z 
4 


m]j = 


ml 


H = 


1 + - 



1 + 



/(A 



2 A2 •"^^ 



2 A2 



(/^l + /02) 



2 A2 



4 A2 



(/(/>1 + f(f>4: + 2 /02) 



1/2 



All the remaining terms in the effective Lagrangian, denoted by dots in eq.( 
terms of the renormalized fields, couplings and masses, as defined in eqs.J 
standard gauge interactions are dictated by the covariant derivative 



(6a) 
(6b) 
(6c) 

(6d) 

), are expressed in 
and eqs.([6]). The 



D, = d, + t^ {T+W; + T-W;)+igz {T' - s^y Q) Z, + teQA,, (7) 

where = (T^ ± i T^)/ 2, Q = T^ + Y, and e = gw sw = 9y cw = Qz cw sw ■ 

Before expressing the HVV interactions of Leg, let us briefly review the observable conse- 
quence of new physics in the gauge boson two point functions in eq.(l5]). First, the ratio of the 



neutral current and the charged current interactions at low energies deviate 



21 



1 



a T 



G 



N.C. 



G 



c.c. 



all ml 
9w/ 



W 



2 A2 ■ 



23| from unity, 
(8) 



Second, the extra kinetic mixing between B^^, and Wj^^^ modifies the 7 and Z boson propagators 



2 A2 



fsw q 



(9) 



in the notation of ref. 22], which contributes to the S parameter 231] 



5 = 16 TT 



s {mz) c {ml 



BW 



(10) 



Here the over-lined couplings e^i^q^), s'^{q'^) = 1 — c^{q^) and g%{q^) are the effective couplings 
that contain the gauge-boson propagator corrections at the momentum transfer 22] • We will 
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examine the constraints on /^^ and fsw from the precision measurements of the weak boson 
properties in the last section of this report. 

The terms describing the HVV couphngs in the effective Lagrangian are now expressed as 

Lf/^ = (1 + c,^^) g m^HW':W-^ + (1 + c,,,) H Z, 



+ 


9z 






mz 


+ 


9z 








+ 


9z 









C2WW HW+ W-^'' + ^1 {d,H) W~ - (d^H) W;^ W+^'' + h.c 



2 ■ 2 

c,z^^H Z^, A^' + c.,z-, ( {d^H) Z, - {d,E) Z^ \ M 









mz - 







2 ^ 



'11^ 



where the 9 dimensionless couphngs, q, parametrize all the non-standard HVV interactions 

,2 



V 



-izz 



^2ZZ 



"277 



■'2WW 



^3ZZ 



4A2 

4A2 



m| 
A2 



~ -^w /bb ~ '^w "^w /bw " "^w Jww 

"^W /bb + 2'-''^ ~ -^w) /bw ~ '^W /wW ) ■^W'^W; 



mz 
A2 

A2 



2A2 



2 2 



/bb + /bw /wW ) "-W ■'W; 

~ /ww ' 



■^W /b /w I ' 



4A2 



2 ^2 



myC 



z^w 



2A2 



/b /w 
" /w 



"^W '^W' 



(12a) 
(12b) 

(12c) 
(12d) 
(12e) 
(12f) 

(12g) 
(12h) 
(12i) 



From the effective Lagrangian of eq. (fTT]) . we obtain the Feynman rule for Vi{pi) — V2{p2) 
H{pfj) vertex as 



THy^^^ipH, Pi, P2) = 9z^z 



f^VlV2 

hi'^' 9t.u + -j-T- P2m Piu 



(13) 



where all three momenta are incoming, pi +P2 + Ph = 0, as shown in the fig. [H Vi and V2 can 
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FIG. 1: The HVV Vertex. 



be {V1V2) = (ZZ), (Z7), (7Z), (77), {W+W-) or {W-W+). The coefficients /^r (Pi, P2) are 

hf^ipi, P2) = 1 + c,zz + ^^^-^2 ^2zz + -f C3^^, (14a) 

lib Til ^ 

hl^iPi,P2) = 2 (c^zz-cszz), (14b) 



for the HZZ couphngs, 



hr{pi,P2) = 2C2,„ (15b) 



for the if77 couphngs, 



^2 I ^2 ^2 ^2 ^2 ^2 



rZ7/ N Pl+P2-^H Pl-P2-^H (.a \ 

„f . d±Ap± _ -p-^r"" C3,„ (16b) 

Z Z 

hpipu P2) = hf{pi, P2) = 2 {c2Zj - Cg^^), (16c) 



for the HZ'-f couphngs. It is to be noted that the HZ'-f coupling c-^^^ has the Feynman rule 
which is not symmetric under an interchange of pi and p2. For the HWW couplings, 



/^r^"(Pi, P2) = /^r"^ (Pi, P2) = hf'^ipu P2) {i = 1, 2) (17a) 
hi [Pi,P2) = [^ + Ciww) COS 9w^ 2 C2ww^ r^sw (l"b) 



^2 {P11P2) = 2 {C2ww-C3ww)- (17c) 

Although we do not consider off-shell Higgs boson contributions in this report, should be 
replaced by {pi +^2)^ in the above Feynman rules when the Higgs-boson is off-shell. 
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III. OPTIMAL OBSERVABLES AND PHASE SPACE 



A. Optimal observables method 

The optimal observables method makes use of all the kinematic distributions which are 
observable in experiments. We therefore summarize our phase-space parameterizations for all 
the Higgs boson production processes in e~^e~ collisions considered in this study, which can be 
generically written as 

e- {h, f ) + e+ (*„ f ) - / (p,4) + / (p- y ) + • (18) 

Here ki and at/ 2 are the four momenta and helicities, respectively, and and Xi/2 are the 




e-ih) f{pi) 

e-(fci) H 

(a) ZH production (b) vector-boson fusion 



FIG. 2: Feynman diagrams for e e"*" ^ f fH. 



four momenta and helicities, respectively, of the produced fermion {i = 1) and anti-fermion 
{i = 2). For / 7^ e, z/g, the processes f llSp occur only through the ZH production diagram as 
shown in figl2](a), whereas for / = e or Ue, both the diagrams figl2](a) and figj2]^b) contribute. 
The effective HVV vertex is depicted by the solid circle in the Feynman diagrams. The ZH 
production process (a) is sensitive to the HZZ and HZ'-f couplings, while the vector-boson 
fusion processes (b) are sensitive to the HWW coupling for f = Ug, and the HZZ, HZ'-f, H'-fj 
couplings for / = e. 

The matrix elements for the processes eq. (fT8l) can in general be expressed as 

n 

M^i^^ = (Msm)^^^^ + (19) 

1=1 
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where {MsM)a\al denotes the SM hehcity amphtude, and q denotes the non-standard couphngs 
of eg. (1121) that contribute to the process. The matrix elements (Mj)^J^^ give the hehcity am- 
phtudes which are proportional to the coupling q. If the e~ and beam polarizations are P 
and P {\P\,\P\ < 1) respectively, the differential cross section can be expressed as 



da{P, P) 



Ssm(P,P;$3) + X^ c^ S,(P,P;<I>3 



i=l 



(20) 



where the non-standard couplings q are assumed to be real and small, and hence the terms 
quadratic in couplings are dropped. Here (i$3 is the 3-body phase space volume of the fjH 

system, and 



(M; 



SsM(P,P;<f3) = ^ 

o'i,o"2,Ai,A2 

gives the differential cross section of the SM. The term proportional to c, 



a,P 



a^P 



SMj 



\A1A2 



(21) 



l+(j^P\ /l + aaP 



2 Re 



(22) 



S,(P,P;$3) = ^ E 

(7l,(T2,Al,A2 

gives the differential distribution which is proportional to q. 

In the optimal observables method, we make full use of the distribution Sj(P, P;$3) in 
order to constrain Cj. For instance, if all Sj(P, P;$3) have different shapes from each other, 
then in principle, we can constrain all the coefficients q simultaneously. For a given integrated 
Luminosity L, the statistical errors of the Cj measurement can be obtained from a function 



N 



k=l 
N 

E 

k=l 



iVi; rP,P)-<(P,P;Q 



NUP^P) 



L Se.p(P, P; $1) A-L [Ssm(P, P; H) + Eti c^ S.(P, P; $|)] A 

v/LSexp(P,P; $^)A 



(23a) 



,(23b) 



where N^-^p{P, P) is the number of events in the k'th bin, and N^^{P, P; q) is the corresponding 
prediction of the theory which depends on the parameters of the SM and Cj. In the second line 
(I23bl) . $3 for A; = 1 to gives the representative phase space point of a bin number k with 
the bin size A. Now, if all the coefficients q are tiny, the experimental result in the k'th bin 
should be approximated by the SM prediction as 



S,.p(P,P; ^1) 
The function can then be expressed as 



^{P,P) 



(24) 



(25) 
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where 



N 



-1 



(26a) 



fc=i 



J:,{P,P; $3) Sj(P,P; $3) 
Ssm(P,P; $3) 



(26b) 



where we take Lq = 100 fb^^ as a nominal integrated luminosity through out this report. If 
the total number of events is sufficiently large, the integral representation in eq. (126bp gives a 
good approximation for the matrix. The value of Xmm ^^"^ ^^e mean value c, depend on 
the actual experimental results, or the small deviation from the equality in eq. (1241) . If the SM 
prediction gives a reasonably good description of the data in most of the phase space region, 
then the statistical errors of q and their correlations are determined solely in terms of the 
covariance matrices V, which is the inverse of the matrix given in eq. (l26|l : 



In practice, however, we should address the following subtleties: 

(i) . If the statistical error becomes small, systematic errors need to be considered. 

(ii) . The results depend on how we split the total Luminosity to different beam polarizations. 

(iii) . Not all the 3-body phase space points are observable in experiments. 

As for the first issue, we assume that the energy and angular resolutions of ILC detectors are 
good enough to justify our integral approximation of eq. fl26|l . and consider only the impacts of 
the luminosity uncertainty as a source of the systematic error which is discussed in section IIXI 
We leave the difficult problem of background contaminations and the spectrum distribution due 
to bremsstrahlung and beamstrahlung photon emissions to future studies. In short, our results 
should be regarded as an ultimate accuracy of the HVV couplings measurement for a perfect 
detector in a background- free environment, when the SM predictions are accurately known. 
On the second point, we provide numerical results for the two very simple cases only: 

1. Unpolarized beam : The total integrated luminosity L is given for collisions with 
(P, P) = (0, 0) at each collider energy ^/s. However, in order to save the length of this 
article, we provide the unpolarized results specifically only for s-channel ZH production 
at y/s = 250 GeV. They are calculated for all the processes at all energy choices and 
are used to evaluate the significance of the beam polarization after all the channels and 
energies are combined in section |Xl 





(27) 
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2. 80% polarized e beam : Exactly half of the total luminosity L is given for collisions with 
(P, p) = (0.8, 0), and the remaining half with (P, P) = (-0.8, 0). 

In general, the covariance matrix depends on the partition of the total luminosity into ex- 
periments with different sets of e~ and e"*" beam polarizations. If the e~ and beams with 
polarizations P = ±|-P| and P = ±|-P|, respectively, are available and if the integrated lumi- 
nosity of L(^pp-j is distributed for each case, then by using the definition of the inverse of the 
covariance matrix for the beam polarization (P, P) in eq. (!26|) . the inverse of the total covariance 
matrix is obtained as 

-^(|P|,|P|) .^-i -^(-|P|,|P|) .^-i -^(|P|,-|P|) .^-i -^(-|p|,-|p|) t^^i /poN 

Our first case is simply L(0, 0) = Ltotai and the second case stands for 

-^(p=o.8,p=o) = -f'(p=-o.8,p=o) = — 2 — ■ (^^-^ 

It should be noted that an equal partition of the total luminosity as above is advantageous for 
the asymmetry measurements, and hence for discriminating among different couplings, while 
L(p=_|P| P=|P|) = Ltotai maximizes the VrVT-fusion cross section. 

We address the third point of limited observable kinematic distributions in the following 
subsection. 



B. Observable phase space for each process 

The observability of the three-body phase space point depends on species of // pair, and also 
on their energy and the scattering angle in the laboratory frame. Let us discuss the following 
4 cases one by one : 

1- f = 

2. f = e 

3. / = Z/e, U^, Ur 

4. f = u, d, s, c, b 

First, for f = ^ and r, all the four momenta pi,p2,PH are measured and used in evaluating 
to estimate the experimental sensitivity. Since only one diagram fig|2]^a) contributes to 
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the cross section, a very tiny fraction of fi^ and escape detection along the beam pipe. We 
therefore use the whole phase space region to evaluate the inverse of the covariance matrix 
eq. fl^B]) . We also assume that the majority of t^t^H events can be made background free by 
selecting the events, in which the r+r~ invariant mass is reconstructed and required to match 
the Z-boson mass. Although this reconstruction is not possible when the Higgs boson decay 
has significant missing momentum, such a.s H ^ and H W^W~, we do not take 

account of resulting reduction of the number of events, since it can be considered as part of 
detection efficiency. 



We do not consider the r polarization in our analysis, because it was found in ref. 13| that 
its impact is not significant once the e~ beam polarization is available. 

Second, for / = e, we can also measure both pi and p2 uniquely. For this process both the 
diagrams figl2](a) and figI2](b) contribute to the cross section because of the possible escape of 
e"*" or e~ (or both) along the beam pipe, the events can be divided into the following four classes 



1. ZH events are selected by requiring 

|cos6'e±| < 0.995, l^e+e- —m^\ <5Tz (30) 

2. Both outgoing are detected (double-tag events) when 

PT(e^) > 1 GeV, Icos^eil < 0.995, \m^+e- -m^] > 5 . (31) 

3. When a photon is exchanged either from electron or positron in the t-channel, the cor- 
responding outgoing tends to escape detection (single-tag events). Those events are 
selected by requiring 

PT{e~) > 1 GeV, I cos ^e- 1 < 0.995 < |cos^e+|, (32a) 



or 



PT(e+) > 1 GeV, |cos^e+| < 0.995 < |cos^e-|, (32b) 



4. When a photon is exchanged from both the legs in the t-channel, both of the outgoing 
tend to escape detection (no-tag event). The selection conditions are 

I cos ^e± I > 0.995. (33) 
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Although the above classification misses a part of the phase space where | cos^e±| < 0.995 and 
PT(e^) < 1 GeV, we find that less than 0.001% (0.02%) of the total e~^e~H events escape from 
this region of phase space at y/s = 500 GeV(250 GeV). 

For the ZH production case (130|) and for the double-tagged eeH events (13T|) . we assume 
that the whole 3-body phase space is observable, and the inverse of the covariance matrix 
is calculated by integrating over the phase space with the selection cuts. In case of no-tagged 
events ( !33|) . most of the events are due to 77 fusion, and it is sensitive to the if 77 coupling. 
We estimate the cross section by using the equivalent real photon approximation in section 
IVIII Since typical transverse momentum of the Higgs boson is smaller than the experimental 
resolution, the only observable kinematic variable is the Higgs boson rapidity 

,„ = iln|^i±^ (34) 

where p^^ is the momentum component along the e~ beam direction. The cross section for 
the no-tag {ee)H events can then be expressed as 

4^ = ^suiVH) + Q T^iiyn) (35) 
dVH Y 

and accordingly, the inverse covariance matrix elements are calculated in terms of ^suiVH) and 
T,i{yH)- In our study, we find in section IVlII that T,i{yH) = constant x SsM(yH), and hence 
only the total production cross section constrains the H'-f-f coupling. 

The single-tag {e)eH events, where either e~ or e"*" escapes detection, are sensitive to both 
H'-f-f and HZ'-/ couplings. We can again use the equivalent real photon approximation for the 
emission from untagged e^, and the differential cross section is expressed as 

da = da{s = sz) D^/e{z) dz , (36) 

where D^/e{z) is the number density of an equivalent real photon of momentum fraction z, and 
the subprocess cross section for '-/e ^ H e can be expressed as 

Ssm(cos^) + ^ Ci Si(cos^) (37) 



dcosO 



where ^ is the scattering angle in the observed e H rest frame. By inserting eq.( l37|) into eq.( l36|) . 
we find 

da = [D^/e{z) Ssm(cos^) + ^^ q D^ie{z) Ei(cos^)] dz dcosO, (38a) 

i 

= [T.sm{z,cos9) + '^^ Ci 'Zi{z, COS 6)] dz dcosO (38b) 
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The inverse covariance matrix is then obtained by replacing S|gM,i}($3) by T,^sM,i}{z,cos9) in 
eq. fl26|) . and by dzd cos 9. 

When / = z/g, z/^ or z/,-, we can measure neither pi nor p2, but only the sum pi + p2 from 
the four-momentum conservation. In this case, only the Higgs boson energy and momentum 
are observable when Higgs boson decays to bb pair, and the observable cross section in the 
laboratory frame is 



da = B{H — i> bb) \T,sM{Ejf, cosOj^) + q Y^i{Ejj, cosdjj)\ dEjj dcosOjj . 



(39) 



(40) 



The covariance matrix is obtained by using the above distributions. In practice, the Higgs 
boson energy from ZH production is peaked in the region 

[s + -2 y/s EjjY^^ - 1712 

and only the cosO^j dependence of the distributions is effective. 

In case of f = u^, in addition, the W boson fusion diagram in figl2](b) contributes. Since 
this amplitude is rather large at high energies, a/s ^ 500 GeV, we make a careful study of 
the py-fusion contribution by using the Higgs-boson rapidity yn and the transverse momentum 
Pth, 



da = B{H bb) 



^suiVH, Pth) + X] ^iiVH, Pth) 



dyu dpTH ■ 



The VT-fusion events are selected by requiring 

{s + mjj-2 y/^ E^Y^^ 



mr 



rrir 



> 5Tz, 



(41) 



(42) 



and in order to avoid contamination with the 77 fusion events, we impose an additional con- 
straint 



Pth > 10 GeV. 



(43) 



When f = q = u,d,s,c,b, we cannot distinguish / from / efficiently. In reference 13j, 
impacts of partial identification of b from b has been studied, and they are found to be negligibly 
small. In this study, therefore, we do not distinguish q from q at all, and we do not distinguish 
quark fiavors. In this limit, the hadronic decays of the Z boson gives two jets with momenta 
Pi and P2 for q and q, respectively, but we cannot distinguish the events from those where g-jet 
has the momentum p2 and g-jet has pi. The appropriate optimal observables are then obtained 
as follows: 



da = - 



Ssm($3) + Ssm($3) + 2^ Ci 



E,($3) + S,(l'a 



d^-i 



(44) 
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FIG. 3: Total cross sections versus e^e~ collision energy ^/s for the five processes that are sensitive to 
the HVV couplings at ILC. All the curves are for rrijj = 120 GeV. The ZH production cross section 
is the sum over all Z ^ ff decay modes with Im^p — mz\ < 5Tz, while the solid thin curves shows 
the Tz = limit. The i^i? H and the double tag e^e~ H events satisfy \'m^j — m^\ > 5Tz, and the 
tagged has | cos^e±| < 0.995 and pTe± > 1 GeV, while pth > 10 GeV is imposed on uuH process. 
The solid thin curves for e^e~ — > VeV(,H and e+e~ — > e~^e~ H give the cross sections calculated from 
the t-channel W and Z boson exchange amplitudes only without imposing the invariant mass cut. 



where $3 is obtained from the phase space point $3 by interchanging pi and p2- The three body 
phase space volume d^s is divided by 2 in order to account for the double counting. The inverse 
covariance matrix is calculated as in equation (1261) . while S's are replaced by S($) + and 
the phase space measure is divided by 2. 
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C. cross section for each process 



Before we start examining the above processes one by one, we present here the total cross 
section of all the five processes for = 120 GeV in figl3] as a function of the cm. energy 
^/s. All the curves show the SM prediction, which are obtained by setting all q to zero. As 
for e~^e~ — > ZH, we plot the cross section for the sum of all the Z boson decay modes when 
\fTijj — mz\ < Sr^. The = limit result is given by the thin curve. The cross section 
for e+e" I'e^eH is obtained by requiring |mj,p — mz\ > 5Tz, eq. fH2l) . and > 10 GeV, 
eq. (j^3|) . in order to suppress the contribution from Z ^ vv and 77 — * H ^ respectively. The 
contribution from the WW fusion (t-channel W exchange) process only without the rHyy cut 
is given by the thin curve, in order to show the relevance of the interference effects and the 
rriyy cut at low energies. The cross section for the process e^e~ — > e^e~H is separated into 
four cases; Z e+e" (|mee — m^l < SF^), double-tag {\mee — > 5F^), single-tag and 
no-tag events, where the final and e~ are tagged when | cos6'e±| < 0.995 and PTe± > 1 GeV. 
The thin curve for the double tag e^e~ H cross section gives the contribution from the ZZ 
fusion only without the rriee cut. The sign of the interference effect is opposite between the 
Vf,V(.B. and e^e~H process, because of the opposite relative sign of the coupling factors between 
the t-channel and s-channel amplitudes. In addition, when i/i ~ "^/z + "^z + 5F^, only the 
rrifj — mz < region contribute with positive (negative) interference in z/gZ/gif [e^e~ H) events, 
while at high y/s the negative (positive) interference from the — ^z > region dominates 
because the magnitude of the t-channel amplitudes grow with Tn^j. 

For the SM cross sections, we use the physical masses, = 80.423 GeV, and = 
91.1876 GeV Q], and the MS couplings at the scale, a(m^) = 1/128.0 and sin^^w("^z) = 
0.2312 for the gauge couplings, except when we use the equivalent real photon distribution 
which is evaluated with a = 1/137. All the calculations are done in the leading order of the 
perturbation theory, since none of our results (errors and correlations) are sensitive to small 
differences in the cross sections. We leave the important task of evaluating the impacts of initial 
state radiation, both from bremsstrahlung and beamstrahlung, for future studies. 

We study the significance of each process in constraining the HVV couplings quantitatively 
at four representative energies, y/s = 250, 350, 500 GeV and 1 TeV. The statistical errors 
are estimated for a nominal integrated luminosity of Lq = 100 fb~^ at each energy. For y/s = 
250 and 350 GeV, we consider only the ZH, Ue^eH and double-tag e^e~H processes, because 
the single-tag and no-tag e^e^H cross sections do not give sufficient number of events in 
the SM. It should be noted that the lyiV-fusion overtakes the ZH production contribution at 
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^/s = 500 GeV, and even the cross section of the double-tag eeH events from ZZ fusion becomes 
comparable to that of the ZH production at ^/s = 1 TeV. The single and no-tag e^e~H events, 
which are sensitive to the HZ'j and H'-ff couplings, respectively, give sufficiently large cross 
sections only at ^/s = 500 GeV and 1 TeV. When we combine results from all the processes 
and from all the energies, we examine the impacts of higher luminosity, by giving L = 500 fb^^ 
at y/s = 500 GeV and at i/i = 1 TeV. 
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FIG. 4: Total cross sections versus for the five processes which are sensitive to the HVV couplings 
at ILC, at (a) ^/s = 250 GeV, (b) 500 GeV and (c) 1 TeV. The tagged has | cos6le±| < 0.995 and 
PTe± > 1 GeV in the laboratory frame. |mjj— < 5 Tz for ZH production and |mjj — m^l > 5 Tz 
for fgPe H and double tag e^e^ H. The thin curves in (a) for vv' H and e'^e~ H show the cross sections 
when the Z ^ ff exclusion cut is removed, and that for ZH shows the Tz = limit, pth 

> 10 GeV 

is imposed on vi'H process. 



The total cross sections are shown in figJH as functions of the Higgs boson mass between 
100 GeV and 200 GeV for (a) = 250 GeV, (b) 500 GeV and (c) 1 TeV. The cross sections 
do not depend strongly on nijj- for < 200 GeV, except for the dominant ZH production 
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cross section at ^/s = 250 GeV which drops sharply from 100 fb at rrifj ~ 150 GeV down to 
1 fb at rrifj = 165 GeV, and vanishes at niH = y/s — mz + 5 ~ 170 GeV because of the final 
state cut \mf f — mz\ < 5 Tz] see eq.( l30l) . The thin curve shows the zero-width limit. The ZH 
production process is dominant up to ~ 165 GeV, above which the VTiy-fusion process 
becomes a main process, and its cross section is above ~1 fb up to ~ 180 GeV. 

In figJUi, we show two curves for the Vei'eH and double-tag eeH processes. The thick 
curves gives the cross sections when the Z ^ ff exclusion cut eg. (142!) is applied, and hence 
they may be regarded as those of WW and ZZ fusion events. On the other hand, the thin 
curves give the total cross sections without the Z ff exclusion cut, and hence they receive 
contribution from the ZH production amplitudes with Z UgU^ or e^e~ transitions. Along 
the thick curves for the WW and ZZ fusion events, the exclusion cut affects the fusion cross 
sections below = 170 GeV ~ ^/s — + 5 F^, which remain almost constant down to 
nifj = 150 GeV ~ y^— — 5 Tz, where the increase of the overall phase space is compensated 
by the increase in the cut-out phase space region. 

At ^/s = 500 GeV in figJlJo, all the cross sections of the t-channel processes increase, and 
that of the s-channel ZH production process decreases. The ZH production and the WW- 
fusion processes have almost the same cross section, over 30 fb up to rrifj = 200 GeV, and the 
cross section of the double-tag eeH event is about one order of magnitude smaller than that 
of Uei^e H . The cross sections of the no-tag and the single-tag eeH events are still small, whose 
maximum is 1 fb and 0.2 fb, respectively, for mjj ~ 2 m^, where the one-loop ^^77 and HZ'^ 
vertices receive the W pair threshold enhancement. 

At ^/s = 1 TeV in figJH^, the cross section of the H^l^-fusion process is above 100 fb, while 
those of the ZH production and the double-tag eeH processes are ~10 fb, almost independent 
of rrij^ up to rrifj^ = 200 GeV. As in figJDD for y/s = 500 GeV, the cross section of the no-tag 
and single-tag eeH processes take the maximum values of 2 fb and 0.2 fb, respectively, around 
rrifj ~ 2 mw- 

IV. e+e" ^ Uei^eH WW-FUSION PROCESS 

The HWW coupling is best measured in WW-fusion production of the Higgs boson, see 
fig|2]3, in the process 

e- [h, ^) + e+ {k,, ^) ^ (^p,, + (^p,, + H{p^), (45) 
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where ki and pi denote four-momenta of each particles, (Ti/2 demote hehcities, and Aj/2 
demote or hehcities. The WW-fusion process contribute only to one helicity amplitude, 
o"! = —a2 = Ai = — A2 = — . Contamination of the HZZ couplings via s-channel ZH produc- 
tion followed by Z ^ uu decays, see figj2^, can be avoided by choosing suitable cuts on the 
observed Higgs boson energy eq. (l4Up . Therefore, the process with "a Higgs + missing energy" 
in e+e~ annihilation can probe HWW couplings independently of the other vector boson cou- 
plings. Once the branching fraction of if ^ WW* decay is known, the measurement of the 
HWW coupling determines the total decay width of the Higgs boson in a model independent 
manner. It then allows us to translate all the Higgs boson branching ratio measurements into 
the measurements of the partial widths, from which we can determine the magnitude of various 
Higgs boson couplings. It is also worth repeating here, that the sensitivity of the HWW cou- 
pling measurement is expected to be better at high energies, because the cross section grows 
as log(s/m^) with energy, in contrast to the ZH production cross section which decreases as 
see figJSl 



A. Helicity Amplitudes and Backgrounds 



The helicity amplitudes for the process given in eg. (1451) can be written as: 



M: 



'I'l = ^-f 5..+ '^A,- 5a.+ rfj^^ Dwit,) Dwih), (46) 



where D]y{ti) = l/{ti — rri^) are the propagator factors with tj = {ki — PiY, (i=l,2) and 



je+i+,+) = ^{P2, +h Y Pl V{k2, 



(47a) 
(47b) 



are the leptonic charged currents. The explicit form of the non-zero components of the massless 
currents in the laboratory frame are given in eqs. (1B1[) in Appendix [Bl The HWW coupling 
can be read off from eqs. (fT3l) and (fTTl) and expressed as 
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where = —p^ (i=l,2). By contracting the t-channel currents with the generahzed HWW 
vertex using eqs. ( 1B2I) in Appendix [HI we can write the hehcity amphtude (H6l) as 

1 



where 



SM 



1 + c 



IWW + ,^2 1 ^2WW 



2 ^ G 

ti+t-2- + — 



4 



SM 



Dw{t2) ^(^1,^2) 
^1 . 



J^{9i,92) = 2 s y/xi X2 cos^ sin ^ , 



^(^1,^2 



X 



2 — xi (1 + cos ^1) + Xi sin 61 cot ^ 
2 — X2 (1 — cos ^2) + ^2 sin ^2 tan — e**^ 



(49) 

(50a) 
(50b) 

(50c) 



Here = 0i — 02, and Xi, ^i, (pi are the energy fraction, polar and azimuthal angles, respectively, 
of z/g (^ = 1) and i/g (^ = 2) in the laboratory frame; see eg. (lAip in Appendix [Al 
The SM contribution to the squared matrix element is 

4-|2 
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SM- 
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Dw'(ti) DH/(t2) ^(^1,^2) 



and the distributions for each non-standard effective couplings are 
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(51) 

(52a) 
(52b) 

(52c) 



As explained in section IIIIl since we cannot observe I'e and Ue momenta, we cannot make 
use of all the distributions Fi in eq. fl52l) to constrain the couplings Ciww- In fact, we can 
measure only the Higgs boson momenta, E^j and cos O^j, or px// and yn- In order to obtain the 
observable weight functions, we perform the integration over the z/g and z/g momenta by using 



the 3-body phase space parametrization in eq.( ]A8cl) as 

,2 



512 TT^ 



d{cose* 



dpTH dvH , 



(53) 



where m^p = s + ■mjj — 2y/s E^, and Ejj = -\/pt/7+""^- "^^^ observable differential cross 
section with initial e~ and beam polarizations P and P, respectively, is 
dV(P, P) 



dyu dpTH 



B{H hh) 



Ssm(P, P; Vh, Pt//) + ^ Q J^cAP^P; Vh, Pth) 



■ (54) 



21 



Since only the left-handed e and right-handed e"*" contribute to the cross section, the weight 
functions are 



(55) 



where i =SM, c^iyiy, Cr^ww ^"^^ ^■mw- "^^^ covariance matrix for the measurements of q 
CiM/vK' ^2ww-- ^zww with an integrated luminosity of L{p^p) is now expressed as, 



^{p,p) 



In the following numerical studies, we set B(if — > 66) = 0.9 for the branching fraction at 
= 120 GeV. 

Before proceeding to the sensitivity analysis of the anomalous couplings, we briefly remind 
ourselves of the potential backgrounds to this measurement. 

1. No tag events in the process e~^e~ — e^e~H via the t-channel 77, 7Z or ZZ fusion 
contribute, but the Higgs-boson cannot have large pth order for the to escape 
detection. We impose a selection cut p-^^ > 10 GeV, eq. (H5]) . which is sufficient to 
suppress the no-tag {ee)H contribution even at ^/s ~ 1 TeV. Contribution from 77 
fusion is negligibly small at ~ a few GeV, as we will show in section IVIIt and that 
from ZZ fusion is estimated to about 5.4%, 4.8% and 3.9% of the signal respectively, in 
the three smallest bins, 10 - 20, 20 - 30 and 30 - 40 GeV at = 1 TeV. 

2. The s-channel ZH production process, where the Z boson decays into a uu pair, can also 
be a background to this measurement at low energies. Since the background events have 
missing mass peaked at m^, we impose the Eh cut, eq. (H2l) . which removes those events 
whose missing mass lies within 5 Tz of the Z boson mass. 

Although small number of background events will survive the above cuts, especially in the 
presence of initial state radiation and finite detector resolution, we expect that their effect can 
be controlled and that they will not affect the main conclusions of the following analysis. 
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B. Sensitivity analysis of the HWW couplings 

1. Sensitivity at ILC-I for ^/s = 500 GeV 

We first perform tlie binned analysis for single and double distribution at -^i = 500 GeV. 
Since the weight functions of eq. (lMl) depend only on yu and Pt^^ calculate the expected 
number of events in a {uh^ Pth) bin with the bin width of AyH = 0.05 and ApxH = 10 GeV. 

We show in figure [5] the weight functions S's integrated over the rapidity, yn- This plot 
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FIG. 5: Histogram showing the pT distributions of the Higgs boson, where the differential cross 
section is integrated over yn in each pth bin of 10 GeV widtli at ^/s = 500 GeV for mjj = 120 GeV. 
^SM gives the SM distribution, and shows the coefficients of the non-standard HWW couplings 




shows that the contribution from the operators corresponding to the couplings, c^ww, Cr^ww 
and C3jyjy can in principle be differentiated by using the distribution, because they have 
different shapes. and has a peak at different Ptw, while S^. changes sign. It 

is not clear, however, how well they can be distinguished. 

The covariance matrix method gives a quantitative answer to this question. Using the p^jj 
distribution, we obtain the matrix elements of the inverse covariance matrix as 

iV-% = BiH ^ bb) Lo ^^(pf\^^j:Pf^ Ap,,, (57) 

^ ^sm{Pth> 

where Pth^ denotes the center of each px// bin. Since the covariance matrix is the same for 
the unpolarized case P = P = 0, and for the polarized e~ beam case with equi-partition of 
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^'^iww — i .15 
^C2ww = ± -16 



the total luminosity into P = \P\ and P = — |P|, as in eq. (l29|) . we show our results without 
specifying the beam polarization. It should be noted, that, because of the simple polarization 
dependence of the WW fusion cross section, as shown in eq. fl55l) . the covariance matrix for the 
general case of asymmetric partition of the total luminosity is simply obtained by adjusting an 
overall normalization factor which is proportional to the total number of events. 
We find for = 500 GeV and Lq = 100 fb"\ 

/ 1 \ 

.9986 1 . (58) 

AcgvKVF = ± .045 V .9989 .9982 1 / 

We observe that the magnitudes of the correlation matrix elements are quite large, implying that 
there is at least one combination of the three couplings which cannot be measured accurately as 
compared to the others. In order to identify which combination of the couplings are measured 
accurately and inaccurately, we obtain the eigenvalues Aj and the associated eigenvector Si of 
the covariance matrix, in terms of which and V are expressed as 

V'-^ = ^ A"-^ ai a-^ , or ^ = • (^9) 

i i 

The resolving power of the measurements can be studied best by showing the eigenvectors and 
their errors ± 

(60a) 
(60b) 
(60c) 

As anticipated, we find that the error of the combination fl60b ) is two orders of magnitude larger 
than those of the other two. In fact, if we sum over the weight functions given in eq. fl55l) . for the 
combination of eq. (l60b ). we find that the sum almost cancel out and hence this combination is 
poorly measured. 

As above, whenever the correlation matrix elements are large, the eigenvectors and their 
errors reveal much more information. Henceforth, we present our results in terms of the eigen- 
vectors and their errors of the covariance matrix V, whenever they are more informative. 

In order to obtain the maximum information from experiments, we further study the two 
dimensional distributions of yn and pth variables. We present the scatter plots of the weight 
functions of Ciww: C2ww and c^wiy in figsl6^, [6)d and [61;, respectively, on the yn, Pth plane. 
We find that adding the rapidity distribution does not help much in resolving the degeneracy. 
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FIG. 6: Scatter plot in the Uh-Pth plane for the coefficient of (a) Cj^^yyy, (b) C2^/^r and (c) c^^^. 
The measure gives {d?'a/dyHdpTu)l^yH^P'TH units of fb for the bin size of pT_f/ = 10 GeV and 
VH = 0.05, at yfs = 500 GeV for m^j = 120 GeV. 



because the yu dependence of the weight functions for c^ry^, C2w\y and c^iyw are all similar. 
Only figjGb gives some hint of an additional independent measurement as it has non-factorisable 
dependencies in terms of yn and Pth, unlike the other two cases. 

By using the double differential distributions, we obtain the optimal covariance matrix from 

U[yH^ Pth ) ^iKVH^ Pth ) 



l, m 



(m)^ 



Ayn ApTjj 



(61) 



where we set the bin size Ay^ = 0.1 and Apx^ = 20 GeV. For the cross section of 77 fb (see 
figure [3]) at ^/s = 500 GeV, we expect 7700 events with Lq = 100 fb~^, when e~ and e"*" beams 
are unpolarized, P = P = 0. The eigenvectors and their errors are 
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When we compare the optimal result (163|1 with the previous one (158|) that used only the p^^j 
distribution, we find more than 30% improvements in the errors of all the couplings and slight 
decrease in the correlations. While by comparing the eigenvectors and their errors of the optimal 
results ( l62l) with (!60|) . we find that the only effect of using the additional information is to reduce 
the error of the least constrained combination by about 30%. Neither the eigenvectors nor the 
errors of the two accurately constrained combinations are affected much by the optimization. 

2. Sensitivity at y/s = 250 GeV and 350 GeV with t-channel only 

The cross section for the WW fusion is small at lower energies; 9.2 fb at ^/s = 250 GeV and 
33.2 fb at 350 GeV as can be seen from the solid thin curve in fig|3l which further reduces to 
5.8 fb and 31.4 fb, respectively, after imposing the Z ^ vv exclusion cut \myy — m^\ > bVz- 
Not only the cross section is small but also contributions from the ZH production amplitudes 
and the interference terms, as well as the effects due to the Z ^ vv exclusion cut, are significant 
at low energies, as can be seen from the dashed thick curve in figOl We first show the results 
of 'theoretical' studies based on the WW fusion events only, which helps us understanding the 
energy dependence of the WW fusion amplitudes when compared with the results at ^fs = 500 
GeV and 1 TeV. A more realistic study will be given in the next subsection. 

At y/s = 250 GeV we will have only 580 events from WW fusion with Lq = 100 fb~^ and 
with 100% efficiency. Thus our integral approximation for the covariance matrix, eg. (156!) . may 
not be reliable. In order to examine the sensitivity of our results on the bin size, we enlarge 
the bin size to Ayn = 0.15 and ^Pth — 30 GeV and count only the contributions from these 
bins with more than 10 events for the nominal luminosity of Lq = 100 fb~^. We find the 
eigenvectors and their la errors at ^/s = 250 GeV with Lq = 100 fb~^ and P = P = to be 

.61 c^ww - -045 C2WW - •'^9 ^^aww = ± -016, (64a) 
•58 q^^y - .65 C2WW + -49 c^ww = ± -029, (64b) 
.54 qyi^^y + .76 C2WW + -^^ c^^ww = ± l-l- (64c) 

We observe that the nature of the eigenvectors are similar to those at -^i = 500 GeV in 
eq. fl62|) . and also the hierarchy between the largest error and the other two persists. The 
overall magnitude of the errors are larger than those of eq. (!62|) at ^/s = 500 GeV by a factor 
of 8, a half of which can be attributed to the reduction of the cross section by a factor of 
13.5 ~ (3.7)^. We confirm that both the eigenvectors and the eigenvalues of eg. (1641) obtained 
with the large bin size do not differ much from those in the integral (infinitesimal bin size) 
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limit. The errors of the two most accurately measured combinations are especially insensitive 
to the details of our binning procedure. 

We also examine the case at ^/s = 350 GeV with Lq = 100 fb~^, since it is natural to study 
the top quark property in detail near the ti threshold. We will have about 3,140 events and 
the eigenvectors and their la errors for P = P = are : 

.44 Ci^y^y - .068 C2WW ~ -89 c^ww = ± -0050, (65a) 
•63 C;^ww ~ -68 C2WW + -37 c^ww = =t -OH' (65b) 
•63 Ci^y^y + .73 C2WW + -26 Cr^ww = =*= -44. (65c) 

Again, we observe the similar hierarchy pattern. 



3. Sensitivity at y/s = 250 GeV and 350 GeV with ZH interference 

Since the cross sections of t-channel WW fusion process are small at ^/s = 250 and 350 GeV, 
contribution from the ZH production amplitudes is rather significant despite the Z ^ vv ex- 
clusion cut \mi,y — > 5Tz- At y/s = 250 GeV, the cross section grows from 5.8 fb to 
6.2 fb mainly because of the 0.7 fb contribution from the s-channel amplitude squared, which 
overcomes the -0.3 fb contribution from the destructive interference between the t-channel and 
s-channel amplitudes. On the other hand, the cross section at ^/s = 350 GeV decreases slightly 
from 31.4 fb to 30.1 fb because of the destructive interference. In addition, ZH production 
followed hj Z ^ t'^t'^ and Z v^-Vj. contributes with 1.4 fb at ^/s = 250 GeV, and 1.0 fb at 
^/s = 350 GeV, after the Z ^ vv exclusion cut. These events are experimentally indistinguish- 
able from the WW fusion z/gZ^e H process. Therefore the total cross sections become 7.6 fb and 
31.1 fb at = 250 GeV and 350 GeV, respectively. 

We show in fig 13 the px// distribution of the SM, S^m, and the weight functions Sj at 
\fs = 250 GeV. The thin lines show the WW fusion contribution only, which are the yn 
integral of the 2-dimensional weight functions used to obtain the results of eq. flM|) and eq. fl65|) . 
The thick dashed lines are obtained after the interference between the t-channel WW fusion 
and the s-channel ZH production amplitudes are taken into account. Since we can safely 
neglect the non-standard contributions in the sub-dominant ZH production amplitudes, we 
can constrain the HWW couplings by using the two dimensional weight functions Sj(?//^,pT_f/) 
of eq. fl55l) . including the interference contributions, where only the SM amplitudes in eq. fl52l) 
are replaced by the sum of the WW fusion and the ZH production amplitudes. It should 
also be noted that the thick dashed SM curve in figH includes the contributions from all the 
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FIG. 7: pt distribution of the weight function for the ui^H process, where the differential distributions 
are integrated over yn at ^/s = 250 GeV for = 120 GeV. SsM gives the SM distribution, and 
Sc, shows the coefficients of the non-standard HWW couphngs Cj = c^\YW'^2WW'^'iWW "^^^ thick 
dashed hues show the weight functions that include both the t- and s-channel contributions, while 
the thin lines are those with the t-channel contribution only. The thick dashed curve for the SM 
distribution, Tism-, is the sum over all the neutrino species (z/g, u^, u^^). 



neutrino flavors, which cannot be distinguished from the signal. The limited phase space at low 
energies gives rise to the complex dependence at 56 GeV and 75 GeV, which are the Higgs 
boson momenta when m^^p = mz + 5 and = rriz — 5 Tz, respectively. We find that there 
are large destructive interference effects for c^y^^Y and Cg^yiy at low pr^. On the contrary the 
^2ww ^3ww curves show strong constructive interferences at high Pth- 

By using the two-dimensional weight functions, we find the eigenvectors and their la errors 
at = 250 GeV with Lq = 100 fb"^ and P = P = to be 

.51 C;^ww + -085 C2WW ~ -86 c^ww = ^ -Ol^, (66a) 
.71 c^ww - -61 C2WW + -36 c^ww = ± -030, (66b) 
■49 Ci^yiy + .79 C2WW + -37 c^ww = ± -51 . (66c) 

Comparing the results of eg. (1661) with those of eq. (l64l) from the t-channel contribution only, we 
find that the errors of the first two eigenvectors are almost the same despite the growth of the 
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SM cross section. It is because the general growth of the magnitude of the SM distribution and 
the weight functions at pt^h > 55 GeV is compensated by the decrease of the weight functions 
at low Pth- The third error decreases because it becomes more difficult to obtain a combination 
of the three weight functions which cancel out. 

Similarly, we present the results for y/s = 350 GeV with Lq = 100 fb~^, where we have 
included the s-channel ZH contribution. 

.39 CivKiy + -032 c^^^ - .92 c^^w = ± -0048, (67a) 
•66 Ci^y^ - .71 c^wiv + -25 c^ww = ± -010' (67b) 
.64 c^yyy^, + .71 c^ww + -29 c^ww = ± -14. (67c) 

Comparing eg. fl67p with eg. (165 p . the errors of the first two eigenvectors are reduced slightly, 
while the largest error is reduced by a factor 3. Again it is a conseguence of the complicated 
interference patterns for the three weight functions. 



4- Sensitivity at y/s = 1 TeV 

Figj3] shows that the WW fusion cross section grows with the cm. energy, and reaches 210 fb 
at ^/s = 1 TeV. Precision studies of Higgs boson properties, including the HWW couplings, 
will be one of the major motivations of the energy upgrade to ILC-II. We find the eigenvectors 
and their errors for Lq = 100 fb~^ and P = P = 0, 

.14 ci^^^ + .038 c^y^y^. - .990 Cg^^^ = ± .00079 (68a) 
.71 Ci^^;^. - .70 C2^^ + .072 cg^^ = ± .0034 (68b) 
•69 c-^^rw + -71 C2WW + -12 c^ww = ± -070 (68c) 

The errors are now a factor of 2 to 3 smaller than those at 500 GeV in eg. (l60l) . The reduction 
factor is bigger than the naive expectation from the cross section ratio (77fb/210fb)i/2 _ 9.6. 
This is because of the s dependence of the weight functions of Cg^y^y and c^ww^ egs. fl52bp 
and fl52cl) . which grows linearly with s. On the other hand, the reduction is not so strong as 
we would expect from the linear growth of the weight functions. We find that the power of the 
weight functions that grow with s for C2^riY and c^^y^y in egs. fl52bp and (I52cp . respectively, is 



greatly reduced because of the non-observability of the z/g and Ue momenta, which results in 
the integration of the factor Q{6i, 62, (p) in eg. (l55l) for the observable weight functions. 

It is instructive to study, just as a reference, how well we could have measured the HWW 
couplings if and z/g momenta were measurable. We would then be able to use the dependence 
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of the weight functions on the full 3-body phase space, xi or X2, cos^i, cos ^2, 4>i — 4>2 in the 
e~^e~ cm. frame; see Appendix |X] for our phase space parameterizations. By using the integral 
approximation of eq. flS^ . we find at y/s = 1 TeV, 

986 = ± -00083 (69a) 

013 c^^,^ = ± .0017 (69b) 
097 = ± -0064. (69c) 

It is remarkable that, no hierarchy among the three eigenvalues survives and that all the three 
couplings are measured accurately and rather independently, as is clear from the single coupling 
dominance of the three weight functions in eq. fl52l) . It is perhaps most remarkable here that 
the error of the most accurately measured combination of the three couplings in a realistic 
environment, 0.00075 in eg. (l68h). is not much different from the corresponding one, 0.00079 in 
eq.dnHK). 



.084 Cmiy^ + .14 C2]yiY 
• 099 c^ww " -986 C2ww - 
.992 c™,y + . 087 C2WW + 



5. Beam Polarization Effects 



All the above results are obtained for unpolarized e"*" and e~ beams, P = P = 0. Since only 
left-handed e~ and right-handed (cri = —0-2 = —) contribute to the WW fusion process, 
the polarization dependence of our results for e~^e~ Vf,Ve can be obtained in a straight 
forward manner. If, for instance, an integrated luminosity of L(^pp-j is devoted to experiments 
with e~ beam with polarization P and e"*" beam with polarization P, then the total covariance 
matrix is simply 



^(p,p) 



;i-p) (i + p) 



^(0,0) 



(70) 



where Lq = 100 fb~^ is the nominal integrated luminosity, adopted throughout this report. It 
is clear from this expression that if we distribute the luminosity equally to all four combination 
of polarizations ; 

T T T T -^total 



^(1^^1,-1^1) 



'{~\p\m 



\~\pi^\p\) 



(71) 



then the total covariance matrix is identical to that of the unpolarized case. The same applies 
for our choice of eq. fl29|) . On the other hand, if we devote the entire luminosity with fixed 
polarization, such as 



'{~\Pl\P\) 



L 



total ; 



(72) 
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then the inverse of the covariance matrix will be larger by a factor of (1 + |P|)(1 + |P|), simply 
because of the large WW fusion cross section. 

V. ZH PRODUCTION PROCESS 

In ZH production via s-channel 7 and Z exchange, we can study HZ'-/ and HZZ vertices. 
Unlike the WW fusion process all the final states are observable when the Z boson decays into 
a pair of charged leptons or quark jets. The technique of optimal observables provides us with 
high discriminating power among the couplings. 

A. Helicity Amplitudes, Backgrounds and Event Selection 

The momentum and helicity assignment in ZH production followed hj Z ^ ff decay is 

{h, y) + (fc2. f ) ^ Z(pz. A) + H(p„), 

U (73) 

We neglect the mass of and the outgoing fermions, and the fermionic chirality conservation 
tells ai = —0-2 = a, and Ai = — A2 = p. Then the helicity amplitudes 

e+e-->ZH-^ffH _ 

JlOL 

where 

jf„ = v{k2, -a/ 2) u{h, a/ 2) (75a) 

jl^p = u{p,, P/2) Y Pp v{p2, -P/2) (75b) 

are the e^e~ annihilation current fl75al) and the Z ^ f f current fl75bp . respectively. Here a 
and /5 denotes the sign of the e~ and / helicity, respectively, and P^ is the chirality projection 
operator P^ = (1 + a 75)/2, sometimes referred to as P_ = Pl and P+ = Pr. The propagators 
and the couplings are denoted as 

D^{s) = l/s, Dz{s) = l/{s-ml + imzTz), (76a) 
9l" = e Qf, g'Jf = gz (rt - Qf 4) , (76b) 
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HZZ 



pI - ml + i niz Tz 



9l'/f2p (74) 



where Q / and Tg^" are the electric charge and the weak isospin of the fermion / with chirahty 
a. 

In the vicinity of the Z boson resonance, fn^j ~ m^, we can factorize the amphtudes of 
eq. fl7i|) into the e'^e~ ZH production part and the Z — > // decay part, summed over the 
decaying Z boson hehcities ; 

M^;-----^/- = E ^'a" D,[pI) T,V^/ (77) 

A 

The production amphtudes T^^ are obtained in the e^e~ coUision cm. frame. Here A = 
(+, 0, — ) is the hehcity of the produced Z boson. By using the generic HZZ and HZ^ 
vertices of eqs. (]14bp and fll6ap . respectively, and by using the e^e~ annihilation currents of 
eq. flB4p in Appendix [B], we find for p\ = m| 



rpZH _ -M 
-'-aX ~ Jla 



e^*(pz,A) (7J 



Ml dlJe^), (79) 



where 



Mr = E 9z V2-S E, Dy{s) (h^ + ^^ 'l , (80a) 

M^ = E 9z V2~s m, Dv{s) /zf^, (80b) 



rA=0 ^ 

y=7,z 



y=7,z 

with (3z = a/I — m| / E"^ and 



h^^-l + r +r i±^+(c -c N s + m^z-m^j, 

z z 
hl^ = -2 [c^zz-c^zz], (81b) 

h!' = 2 C3,, + {c,z, - C3Z7) '^"V"^ ^ (81c) 



hp = -2{c,z,-c,z,)- (81d) 
In eq.(179l). the d functions 

dl.-.M--^, rfU.W = ii^ (82) 

dictate the overlap of the initial (e^e~) and the final (ZH) state angular momentum states 
in the J = 1 channel, in terms of the opening angle 6 between the electron and the Z-boson 
momenta in the e~^e~ cm. frame. 
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The Z boson decay amplitude is expressed in the Z ^ f rest frame simply as, 

T!f" = gl" eP{pz, A) = ^ m, 4,^(r) (83) 

for Ai = — A2 = (3 massless quarks and leptons, and again in the zero-width limit, p| = m|. 
Here 9* and 0* are the polar and azimuthal angles of the fermion / in the decaying Z boson 
rest frame, where the polar axis is chosen along the Z boson momentum direction in the e^e~ 
rest frame, and (p* is measured from the e+e~ ZH scattering plane; see Appendix \K\ 

From the above helicity amplitudes, we can calculate the weight functions Sc^. The dif- 
ferential cross section, with the e~ polarization a, summed over the final state helicities (3 is 



1 






/3 A 

We can simplify the phase space (i$3 in the small width limit as 



TT 5{pl-ml) dpi 
rrizTz 2 vr 



2 7r 

d<^2{Z^ff) d^2{ZH) 



2mzTz 

Here (i$3 is the reduced 3-body phase space element, 

d$3 = d^2{Z^ff) d<^2{ZH) = {Sny /^(^' ^1 dcosOz dcosO* 

where P{a,b) is given by eq. (lA9aP in Appendix Rl 
The differential cross section is now expressed as 



A, A' 



B4a) 
^4b) 



(85a) 
(85b) 
(85c) 



(86) 



In this equation, we introduced the Z ^ f f decay density matrix by summing 
helicities and colors of / and /, 

1 



(87) 
over the 



Piy 



rpZ^ff IrpZ^ff 



■A'/3 



color P 
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which is normahzed as 



AA' 



- (89) 
8 vr Tz 8 TT ^ ^ 

By using eqs.f lB5P - flB7p in Appendix [B] for the Z ^ ff currents and the Z boson polarization 

vectors, we obtain a compact expression for the Z boson decay density matrix 



Pix' 



12 TT B 



/ 



v 



i+cos^ e* 


sin d* cos d* 


2 






sin 9* cos 




sin^ e* 

sin 6* cosS* ^— 


V2 
sin^ d* - 


-2i<t>* 


2 ^ 








^ — COS 


n* sinO* i 


+A/ 


sin 61*^ 






\ 


sin 6** — 
V2 ^ 



i 0* 



sin^ 8* 2i< 
2 ^ 
sin 6** COS0* ^ 

1+cos^ e* 

2 





sin 6* (f)* 
V2 ^ 



(90) 



for A, A' = 0, 



B 



Here Bj is the Z ^ ff decay branching fraction 
Viz ^ //) _ Nf 



Tz 16 TT Tz 

with Nf = 1 for leptons and Nf = 3 for quarks, and 

(/J'Y - id") 



{A"f + {<,l"f 



A, 



(91) 



(92) 



is the left-right asymmetry parameter of the Zff couplings. In the following analysis, we use 
the tree-level expression for the asymmetry parameter eq. (1921) for the charged lepton decays, 
whereas we replace the branching fraction eq. (l9T]) by the observed values 21] in simulating the 
number of events. 

Let us now obtain the weight functions for polarized beams. The differential cross section 
for polarized e"*" e~ beams is expressed as 



daf (P, P) 



(93) 



where the weight functions are obtained by summing over the e helicity a, 

Q, \ / \ / 

E (^sm)^ (MsM)f dl,,i9,) p{. 



(94a) 



A,A' 



-y. 

2 S ^ 



1 + a P 



A,A' 



1-a P 
2 

, A 



{MX (MIm): + (Msm): {M*X <a(^z) diAOz) piy 



. (94b) 
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Here the ZH production matrix elements, (Msm)^ and (M^-)^ = (Mcizv)^' obtained from 
eqs.(l80l)-(l8Tl): 



v2 s 



s — m 



— 9z 9a 7a, 



{Mc,,X = iMsM)l Rv{s) fUs) {k = 2, 3), 



(95a) 
(95b) 



with 

7± 
fUs) 
where 6n 



7o 



nir 



r: 



cot 9 



w 



s — m 



nir 



s — 

s 

2 



" — cIao 



mr 



7o 



, . 2 7^ - 1 

— ^ + 0), 



J AO 



■^00 — ^77 

expressed as 



z ^z 7o 

= 5zz = 1 and 6±o = 6z-y = S^z 







-- tan 2 6'w 




-m| 




-R^ - 
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/2^ 
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» = 


s + 


1 - 


f 5ao ■ 



2 72 - 1 



7o 



2v^ 7o' 



(96a) 
(96b) 
,(96c) 
(96d) 



j:Up, p- $3) 



^\ 2 ) \ 2 ) {s-mlf 



^z 7o 

0. The weiffht functions T^l are now 

^ee\2 ^2 



(97a) 



T.l (P, P: 



A, A' 



E 



l+aP\ /l-a/>\ 9l(gf") 



2 y V 2 y (s-m|)2 
E 7A' (/.V(^) + fkUs)] dUOz) dlAOz) Mpix') 



A, A' 



(97b) 



for A; = 2, 3 and V = Z,'y. The weight function for c^zz simply twice the SM one, eq.(l9j 
Since all our matrix elements are real, only the real parts of the density matrix elements pj^y 
in eq. (l90l) contribute to the weight functions. 

The weight functions depend both on the beam polarization and on the final fermion species. 
For brevity, we give expressions for covariance matrices for electron polarization P, with un- 
polarized positron (P = 0). For / = e, /i or r, we can use the full phase space information of 
6, 9* and 0*, and accordingly the covariance matrix is obtained as 



S^P,0; l>3) Sl (P,0; $3) 
' S^^(P, 0; $3) 



(9^ 
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where Lq = 100 fb'^ 

On the other hand, when the outgoing fermion is invisible, i.e. , for / = z/g, t'^ or Ur, we can 
observe only the cos 9^ = — cos 9z distribution, and the covariance matrix is. 



(Vp) = B(H ^ bb) Ln / dcos9r 

^ ' ' jT.l^{P,Q- ^^)d cos 9* d(t)* 



(99) 



Here again, we multiply the branching fraction B(iJ bb) = 0.9 because the H ^ tt decay 
mode may not be useful in the analysis. 

When the outgoing fermions are quarks, f = u,d, c, s, b, we use only the half phase space 
information for 9* and 0* because the g-jet and g-jet cannot be distinguished. We find 



i:i{P, 0; $3) 0; $3)) (s^'^(P, 0; <l>3) + (P, 0; $3)) . 

= Lo I ^ y ^ -d^3 

's^sm(^, 0; <I3)+S^sm(^, 0; $3) 



2 

(100) 

where $3 is the phase space obtained by interchanging g-jet and g-jet; see eq.f l44p . Finally, the 
electron beam polarization is taken into account as 

P=±0.8 " f=l,t^,q ^ ^ 

where Lp is the integrated luminosity for each electron beam polarization. 

We expect no significant background to this channel. But the WW fusion process, e^e" — > 
Vf, Vf. if, can mimic the HZ production when the Z decays into neutrinos. This background is 
avoided by demanding the missing mass to lie within five times the Z boson width, see eq. fHU]) . 
After the cut, we confirm numerically that our results are not affected by the interference 
effects, in the UgUeH channel. 



B. Sensitivity analysis of HZZ and HZj coupling 

We first note that the five CP conserving dimensionless effective couplings c^^^, Cg^^, c-^^^, 
C2z^ and c^z-i coming from the dimension-six operators, contribute to the process only through 
the four form factors of eq. flSlj) . In particular, only two combinations of the three couplings 
c^zzi ^2zz ^"^^ ^3zz appear in the from factors of the HZZ and HZ'-j vertices. Consequently, 
we cannot determine the couplings independently at a fixed collision energy. 
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In the following analysis, we show the eigenvectors for the four combinations of the effective 
couplings, 



^2zz — ^2zz "I ; 2 ^izz^ ^■izz — ^2,zz "I ; 2" ^izz^ ^22^, and Cg^^, (102) 
s + s + ' ' 

which can be measured simultaneously at each energy. Note that the unmeasurable combination 

Tn?' 

^\ZZ ~ {(^2ZZ + Cszz) I 2 (103) 

changes with the collision energy ^/s, and hence moderate constraint on the coupling c-^zzi 
which measures the strength of the standard HZZ coupling, can be obtained after combining 
experiments at different energies as discussed in ref. 13 1. 

The covariance matrices are expressed in terms of the eigenvectors and the errors. With 
un-polarized electron beam, \P\ = 0, we find for -^i = 250 GeV and L = 100 fb~^, 

.88 c'^zz + -41 c'^zz + -097 c^^^ + .15 Cg^^ = ±.00048, (104a) 

.45 d^zz - -86 d^zz - -026 c^^^ - .24 c^^^ = ±.0042, (104b) 

.068 d^zz + -28 d^zz - -54 c^z, ' -79 c,z, = ±-0054, (104c) 

.046 d,zz + -10 d,zz + -83 c,z, - -54 c,z, = ±.072. (104d) 

The smallness of the magnitudes of the coefficients of the HZ'-f couplings in the first two 
eigenvectors and that of the coefficients of the HZZ couplings in the latter eigenvectors reveals 
that the HZZ couplings and HZ'-f couplings are measured rather independently. Inspection 



of the first two eigen- vectors in eg. (11041) . which are most tightly constrained, tells that the 
couplings d2zz ^i-nd d^zz '^^^ t)e measured much more accurately than the other two couplings, 
C2z^ and Cg^^, and also their sum in eg. (1104b .) can be measured almost one order of magnitude 
better than their difference in eg. (1104b ). The errors of the HZ'-f couplings, C2Zj and c^z-y^ which 
can be deciphered from the third and the fourth eigenvectors, are significantly larger; those of 
their sum in eg. (1104b ) and their difference in eg. (ll04d ) are, respectively, more than one and 
two orders of magnitude larger than the smallest error in eg. fl 104b ). This can be attributed 
to the cancellation of the contributions from the HZ'j couplings between the left handed and 
right handed electron contributions, where the interference term between the 7 and Z exchange 
amplitudes change sign, 

gr ^ + gr gn' = - e {gl^^ + g'n') = e gz (1/2 - 2 sin^ Ow) <^ e gz . (105) 



Expressing eg. (11041) in terms of the covariance matrix, we reproduce identical results given 
in eg. (5.4) of ref. [131], even though we do not take into account r-polarization and partial 
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identification of botli b and b jets. Small improvements of the covariance matrix due to these 
additional measurements can be hidden in the numerical uncertainty of the correlation matrix 
elements in the standard expression in terms of the errors and their correlations. 

The constraint on the couplings become significantly stronger when we introduce the e~ 
beam polarization. With |P|=80%, L(p=o.8) = -^(p=-o.8) = 50 fb~^ at ^/s = 250 GeV, we find 



.056 c'^zz + -027 d^zz + -53 c^z^ + .85 c^z^ = ±.00024, (106a) 

•90 c'^zz + -43 c'szz - -032 c^z^ - .053 Cg^^ = ±.00049, (106b) 

-.094 c'^zz + -20 c',zz " -83 c,z^ + .51 c,z^ = ±.0034, (106c) 

.42 c',zz - -88 c',zz - -19 c,z^ + -12 c,z^ = ±.0043. (106d) 



It is salient that the errors of the couplings Cg^^ and c'^zz^ that of their sum in eq. fll06b ) and 
that of their difference fll06d ). do not improve at all by introducing the beam polarization, 
from the corresponding ones in eqs. fll04b ) and fll04b ). On the other hand, the errors of the 
two eigenvectors with dominant Cj^^ and c^z^y components are reduced by a factor of 20. This 
observation can be explained as follows: 

1. The C2ZZ ^zzz coupling contribution depends on the beam polarization exactly the 
same way as the SM. In particular, the weight functions for C2zz ^^^d '^zzz eq. (197b[) 
have the same polarization dependence as the SM contribution in eq. fl97ap . essentially 
because = 1 in eq. (]97bp . 

2. When a photon propagates in the s-channel, the c^z-^f and Cg^^ contributions to the 
amplitude do not change by the beam polarization since the electromagnetic interactions 
are chirality blind. But the interference of the 7 exchange and the SM Z boson exchange 
amplitude changes sign; see eq. fll05p . 

We present our results for = 350 GeV for L = 100 fb"^ and \P\ = 80%. W 



.029 C2ZZ + -030 c'szz + -34 c^z^ + .94 c^z-y = ±-00015, (107a) 

.71 c'^zz + -70 c',zz - -013 c^z^ - .040 c,z, = ±.00040, (107b) 

.67 c'^zz - -68 c',zz - -28 c^z, + -100 c,z^ = ±.0014, (107c) 

-.21 c'^zz + -21 c.,zz - -90 c^z^ + .32 c^z. = ±-0012. (107d) 



[*] From now on, wc do not report our results for |P| = in order to save space. Individual results with |P| = 
will be available on request from the authors. 
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We first note that the eigenvectors and their ordering remain similar to those at ^/s = 250 GeV. 
We observe that c'^zz more accurately measured at high energies where the amplitude for 
the longitudinal Z boson (A = 0) grows, mainly because of the 70 = Ez/rriz factor in eq.(l95b). 
The weight functions for the couplings d^zz ^^'^ ^zz-y grows as 7q from the the dominant A = 
amplitudes; see eq. (l96l i). In contrast, the weight functions for c^zz Cg^^ do not increase 
at high energies, since the A = amplitude tends to be canceled by the A = ± (transverse) 
modes; see eg. (196b). This growth of the weight functions at high energies explains the general 
decrease of all the errors at -y/s = 350 GeV in eg. (11071) as compared to those at y/s = 250 GeV 
in eg. (11061) . despite the decrease of the total cross section from 260 fb to 130 fb. 
At ^/s = 500 GeV for Lq = 100 fb"^ with P = ±0.8, we find 



.015 C2ZZ 


+ -035 c'^zz + -18 C2z^ 


+ .98 


^3Z7 


= ±.000099, 


(108a) 


38 C2ZZ + 


.92 c'.zz - -0085 C2z^ - 


- .037 


C3Z7 


= ±.00028, 


(108b) 


.24 C2ZZ 


- -097 c'szz + -95 C2z^ 


- .17 




= ±.00086, 


(108c) 


.89 C2ZZ 


— .37 o^zz " -25 c^z-f - 


- .046 


C3Z7 


= ±.0010. 


(108d) 



Despite smaller cross-section at ^/s = 500 GeV, all the errors in eg. (11081) are smaller than those 
in eg.([I07D for = 350 GeV. 

The energy dependence of the measurement at higher energies {^/s ^ rriz) can be read off 
from Table [H where we show the high energy limit of the matrix elements (McJ^ of eg.( !95l) . 
Since the weight functions are the interference between the SM amplitudes and the Cj ampli- 





Zt(\ = ±) 

Energy Dependence coefficient 


Zl{\ = 0) 

Energy Dependence coefficient 


SM 


1/70 gz at' 


1 9z at' 


^2ZZ 


70 4 gz gt' 


1/70 az at' (5- "4/"^!) 


'^3ZZ 


1/70 azat'^ml/ml 


70 4 gz gf 


^2Z-i 


70 4 e 


1/ 70 e (5 - m]j/ m|) 


'^3Z7 


70 4 e 


70 8 e 



TABLE I: The high energy hmit of the matrix elements (Mc.)^ in eq. (|95p . Here 70 = Ez/ n^z- 



tudes, we immediately find that c^zz ^"^^ ^zz-y measured most accurately at high energies 
from the Zl contribution, because their weight functions grow linearly with 70 = Ez/mz- The 
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couplings C2ZZ ^^^d c^z^ are best measured from the contribution, whose weight functions 
do not vanish at high energies. 



(109a) 
(109b) 
(109c) 
(109d) 

All the four couplings are now measured rather independently, as we can observe from the single 
coupling dominance in all the four eigenvectors. The error of Cg^^ and d^zz becomes half of 
those at 500 GeV in eq. fllOSp . because their weight functions are proportional to 70. On the 
other hand, the errors of d2zz and Cg^^ do not improve because the weight functions remain 
constant at high energies. 



picbcllL Lilt; IcbUlLb lUl 'y — 


1 TeV for Lq = 100 fb~ 


-1 ^xnih P — 
WlLll 1 — 


.004 + .039 4^^ + 


.047c2^^ + 


.998 C3^, = 


+.000046 


.092 c^zz + -995 c^zz ~ 


■003 C2Z,- 


■039 c,,^ = 


+.00013 


-.16 + .014 4^^- 


.986 C2z^ + 


M7c,z, = 


+.00086 


.98 C2ZZ ~ -092 Cg^^ - 


- -16 022^ + 


■007 c,,^ = 


+.0012. 



VI. e+e- e+e-H DOUBLE-TAG ZZ-FUSION PROCESSES 

In this section we study the HZZ, HZ'-f and H'-f-f couplings via t-channel Vector Boson 
Fusion (VBF) process, e~^e~ e'^e^H; see fig|2b. The merit of the t-channel VBF processes 
is that the cross-section grows with ^/s, and also, we expect sensitivity to the coupling c^zz^ 
which cannot be measured independently from the other couplings in the s-channel process, 
e+e- ZH. 



A. Helicity Amplitudes and Background 

The momentum and helicity assignments are the same as those in eg. (145!) for the W^-boson 
fusion process : 

e- {h, ^) + e- {k2, f ) ^ e- (p„ ^) + e- (p^, ^) + Hip,). (110) 

Here again we neglect the mass of e^, and only two diagrams of fig|2|i {ZH production) 
and figOD (VBF) contribute. The electron chirality conservation tells = A,/ for the 
VBF amplitudes, and ai = —02 and Ai = — A2 for ZH production amplitudes. When 
0\ = \\ = —02 = — A2, the two amplitudes interfere. The ZH production contribution is 
suppressed by requiring the invariant mass of the e+e~ pair to be away from m^, eq. fH21) . 
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The helicity amplitudes for the VBF process are given as 



J 



e (cti,Ai) 



■yee 
-a- 2 



■yee 

-(T2 



(111) 



for Ai = (Ji, and A2 = o"2- The propagator factors Dv{U) and the Vee couphngs g^'^'^ are given 
in eq. fl76|) . The t-channel currents 



Je-iaa) = ^(^1' "/ 2) ^ Pa U^k^ «/ 2) , 

Jlnm = ^(^2, P/2) Y P-p vip2, (3/2) 



are evaluated in the laboratory frame, 



Je (aa) 



Je+m 



Js xi COS—, sm— e -% a sm — e cos — 
V i I 2 ' 2 2 2 



( sin^, cos^ e-^*^^ 2/5 cos ^ e~'^"^^ 



02 

2 



2 



0, 



sm 



(112a) 
(112b) 

(113a) 
(113b) 



2 ' ' 2 ' 2 

in the rrie = limit, where the four- momenta of the final electron and positron are parametrized 
as in eq. (lAip . The helicity amplitudes are then expressed as 



m: 



a(3 



af3 



gzm, E gI'''P>vAti)P>v2it2)g':TJ'{0i,92) 

Vi=Z,7 V2=Z,-y 



rrir. 



, (114) 



where the functions JF, 7Y, Q and Q' are given in eq. (]B2I) and ( 1B3I) of Appendix [Bl The above 
helicity amplitudes can be expressed as 



(115) 



where 



(M, 



(Msm)"? 

^?Z^Z^7f"^7-F^^(^l)^z(t2) ^(^1,^2) I 5a;3 m,^2,0) ),(116a) 



iMc2.,V2)f, = (^Sm)S 



{Mi 



a/3 -Dyi(ti) -0^2(^2) 



a/3 -Dyi(tl) -0^2(^2) 



/2V1V2 + 
/3V1V2 - 



g (^a/j Q' + ^a -/3 ^ 

4m| 5^^ n{9^,92,4>) +Sa,-f3 

S Sgp G' + ^a,-/3 g 

4m| 6ap n{9i,e2,(p) +Sa,-p 



(116b) 
(116c) 
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with 



fizz — 1, /2ZZ — /2Z7 — /27Z — /277 — — J ^, (117a) 



fzzz — — f-, /3Z7 — — T fs-fZ — — T -■ (117b) 



mjj mjj -ti+t2 mjj + ti - ^2 



rriz rriz 

It should be noted that in the notation of eq. flllSp . the coefficients of the HZ'^ couphngs should 
be interpreted as 



CkZ^{M^^^X^p + Cfc7z(Me,^^)°g = CfcZ7 (^c.^jaj^ + [M^^^^fJl^ 



:ii8) 



for /c = 2, 3, by using the expressions in eq. flll6b ) and flll6b ). 

The SM distribution Ssm and the weight functions Ec- depend on the 3-body phase space 
in the laboratory frame which is parametrized as in eq. flASep . 

^^3 = 77^;^ r ' \ ' ^^YTon dcf), (119) 

1024 TT* [1 — xi (1 — cos 6^12)/ 2J^ 

where = 02 ~ 0i- The differential cross section with e~ and e"*" beam polarizations, P and P, 
respectively, is expressed as in eq. (l20l) where the weight functions are 

2.m(p. p; 4.) = ^ E (^) (^) l(^^sM)Sl^ 

EJP.P; *3) = ^ E (iif^) (ii^) 2 Re [(M,):; (MsVOS] ,(120a) 

The six non-standard couplings, q = c^^z^ '^2ZZi '^zzzi ^2^75 ^3^^ and C2^^, contribute to the 
t-channel Z and 7 exchange processes, and the inverse of the covariant matrix is then evaluated 
as in eq. fl26bl) . As shown in figl3]and noted in section [III C[ the interference contribution from 
the s-channel ZH production amplitudes is not negligible at \fs = 250 GeV and at 350 GeV. 
We therefore replace the SM amplitude in eq. fll20p by the sum of the t-channel and s-channel 
amplitudes at these energies. That is, we neglect contributions from the anomalies in the 
suppressed ZH production amplitudes. 



B. Constraint on the HZZ, HZ^ and H^^ couplings 

1. Sensitivity at y/s = 250 GeV and 350 GeV with the integration approximation 

We ffist present the 'theoretical' result for the contribution from the t-channel vector boson 
fusion, using the standard integration procedure for ^/s = 250 and 350 GeV. Later on we take 
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into consideration the interference effects from s-channel ZH production and adopt a more 
realistic approach to evaluate the inverse of the respective covariance matrices. 

Our results for the purely t-channel process at i/i = 250 GeV for Lq = 100 fb~^ with 
\P\ = 80% are : 



001ci^^-.000c2^^-. 


^O'^c^zz - 


.001 C2z,+ 


.99997 Cg^^ - 


.008 C2^^ = 


±.0014, 


(121a) 


.028 c^zz - -OQ^c^zz - 


V. 00b c^zz 


- .065 C2^^ 


-.008C3^,- 


.995 C2^^ = 


± .0040, 


(121b) 


.009 c,zz- -097 c,zz - 


KOlOCg^^ 


+ .993c2^^ 


+ .000c3^^- 


.058 C2^^ = 


±.0099, 


(121c) 




— .79 C322 


+ .056 c^z^ 


-. 002 c^z-y - 


.036 C2^^ = 


± .047, 


(121d) 


.73 Cizz ~ -65 C2ZZ 


— .19 c^zz 


-.OUc^z, 


-. 000 c^z,+ 


.068 C2^^ = 


±.067, 


(121e) 


.62 Cizz + -52 C2ZZ 


+ .58 c^zz 


+ .039c2^^ 


+ . 000 c^z,- 


•017c2^^ = 


±.17. 


(121f) 



The cross section for the t-channel double-tag eeH process is only 0.83 fb as can be seen from 
the solid thin curve in figJH Imposing the exclusion cut |me+e- — > bVz, it reduces to 
0.55 fb, or 55 events for 100 fb~^. Nevertheless, we present the above results in the integral 
approximation (the small bin size limit) over the 3-body phase space, as a reference to study the 
energy dependence. Binning effects for small statistics is reported in the next sub-subsection. 

It is remarkable that the two HZ'-) couplings and the if77 coupling are rather uniquely 
constrained in the first three eigenvectors of eq. fll2ip with small errors. This is a consequence 
of the combined effect of the strong e~ beam polarization dependence of the Z — 7 interference 
term (even though it is effective only for the 7-exchange between e~ and H) and also from 
the kinematic difference between the 7 and Z propagator factors when or \t2\ or both are 
significantly smaller than m|. In fact, we confirm the latter effects from the unpolarized beam 



(P = 0) case where the three couplings are still measured rather uniquely [4^; C2^^ is now 
constrained most accurately with the error ±.0041, c^z-^ has an error of ±.011, while C2^^ is 
poorly constrained. The three HZZ couplings are constrained rather weakly in the last three 
lines of eq.ffml). 
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Similarly at = 350 GeV for L = 100 fb"^ with \P\ = 80%, we find 



002 c^zz - -000 C2ZZ - -003 c^zz + -001 c^z.^ + -99998 Cg^^ - 


-.006c2^^ = 


± 


.00085, 


(122a) 








.0035, 


(122b) 


.025 c^zz - -099 c^zz + -006 c^zz + -98 c^z^ - .003 Cg^^ 


— .19 C2^^ = 


± 


.0053, 


(122c) 


.29 c^zz + -41 - -86 c^zz + -030 c^^^ - .004 c^z^ - 


- -04702^^ = 


± 


.016, 


(122d) 


.47q^2 - -82 02^^ - .24 Cg^^ - .057 Cg^^ - .000 Cg^^ 


+ -19c2^^ = 


± 


.025, 


(122e) 


.83 c^zz + -34 c^zz + -44 c^zz + -004 c^z^ - .000 Cg^^ - 


- .032 02^^ = 


± 


.054. 


(122f) 



Here the cross section is 2.9 fb which is almost a factor 3.5 times larger than those at 250 GeV; 
see figlSl The effect of the exclusion cut is negligible and the cross section is reduced to 
2.8 fb which would mean 280 events with 100 fb^^. Therefore, we naively expect a factor of 2 
improvements in the statistical error when compared with the results eq.f ll2ip . The results in 
eq. (1122p show that the error of the HZZ couplings decreases to about 1/3, those of the HZ'y 
couplings to ~ 1/2, while the couping does not show a significant improvement. The if77 
coupling measurement does not improve because the amplitudes with < m| do not increase 
much with ^/s. We notice that the two couplings c^^^ and Cg^^ now have some correlation, as 
can be seen from the combinations (0.21c22^ + 0.96c2^^) in eq. (]122bp and (0.98c2^^ — 0.19c2^^) 
in eq.f ll22c l ) . It is because of their similar behavior when both and |t2| are larger than m^. 
Such region of the phase space is tiny at ^/s = 250 GeV, but starts appearing at 350 GeV, and 
will be dominating at higher energies. 

2. Sensitivity at y/s = 250 GeV and 350 GeV with the table method including s-channel contri- 
butions 

As the t-channel cross sections are very small at ^/s = 250 GeV and 350 GeV, the interference 
effects due to the contribution from ZH amplitudes along with the exclusion |me+e- — mz\ > 
5 Tz can be significant as can be seen from the thick dashed line, in figlSl To simulate the 
realistic experimental situation we adopt the table method to calculate V~^, which is more 
realistic when the expected number of events is not large. Here we examine the two dimensional 
weight functions in terms of the momentum transfers |ti| and |t2|- 

At ^/s = 250 GeV with Lq = 100 fb~^ and \P\ = 80%, the cross section of the double tag 
e~^e~ H process becomes 0.93 fb, after including the contribution from the s-channel amplitudes; 
see figj3l Because of the smallness of the signal events (93 events for 100 fb^^ with 100% 
efficiency), we integrate the weight functions over the azimuthal angle and rrie+e-, and divide 
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|ti|'s into 3 regions each as follows: \ti\ < (45 GeV)^, (45 GeV)^ < \ti\ < (90 GeV)^ and 
\ti\ > (90 GeV)2. With this binning, all the 18 bins (9 bins each for P = 0.8 and P = -0.8) 
have approximately 5 events. We obtain the matrix by summing over the contributions 
from 18 bins, and find 



.006 c^zz 


•0 '^2ZZ 


-.010C3^^ 


+ -31 C2z^, 


+ .95c3^^ + 


.051 C2^^ = 


± .0043, 


(123a) 


moc.zz 


+ -13 C2ZZ 




+ -93 c^z^ 


- -31 C3^7 + 


.065 C2^^ = 


±.013, 


(123b) 


.2Ac^zz 


+ -63 C2ZZ 


+ -60 c^zz 




+ .064c3^^- 


- .41 C2^^ = 


±.016, 


(123c) 


.2Qc^zz 


+ . 29 C2ZZ 


+ -18^32^ 




-.004c3^^4 


- .89 C2^^ = 


± .043, 


(123d) 


.80 c^zz^ 


' -13 C2ZZ ~ 


- .56 c^zz ^ 


-.010C2^^ 


-.005C3^,- 


- .16c2^^ = 


±.088, 


(123e) 


48c,zz- 


.70 c^zz + 


.53 c^zz + 


.012 C2^^ - 


-.001C3^,- 


.015 C2^^ = 


±.37. 


(123f) 



In spite of the increase of the cross section, from 0.55 fb to 0.93 fb, we observe that the errors of 
all the eigenvectors in eg. (11231) are larger than these in eg. (11211) which are obtained by using the 
integral approximation. Especially, the error of the eigenvector with dominant c^^^ component 
in eg. fll23dp becomes an order of magnitude larger than the corresponding one in eg. fll21b|) . 
The main reason for the enhancement of the errors is the loss of information due to the large 
bin size. For instance, the error of c^^^ is reduced significantly if we divide the < (45 GeV)^ 
bin into 2 bins, while that of 03^^ is reduced if we divide the > (90 GeV)^ bin into 2 bins. In 
the very high luminosity limit, we can decrease the bin size, and the integral limit of eg. fll2ip . 
corrected for the interference effects, is obtained. In the following analysis, however, we will 
use the result eg. (11231) as the contribution from the double-tag eeH events at = 250 GeV 
for 100 fb~^ with \P\ = 0.8. 

We repeat the same exercise at ^/s = 350 GeV for Lq = 100 fb~^ with \P\ = 0.8. As can be 
seen from figj3]the cross section increases to 3.3 fb from 2.9 fb due to constructive interference. 
Here we considered 6 bins each for and |t2|, to make the number of event almost the same 
in all the 72 bins. We find : 



004ci^^± 


.005 C2ZZ- 


•002c3^^± 


.llC2z^ + 


.99 c^z^ + • 


014 C2^^ = 


±.0012, 


(124a) 


.007 c,zz^ 


- .006 C2ZZ - 


-.072c3^^- 


^•97c2^^- 


- •IIC327 + 


•22 C2,, = 


± .0074, 


(124b) 


.21 Cizz 


± .31 C2ZZ ^ 


-•44c3^^± 


•21 C2z^i — 


•013 C3^^ - 


• 79 ^2^^ = 


±.011, 


(124c) 


.67 c^zz ^ 


- .20 C2ZZ " 


.70 c^zz ~ ■ 


02QC2Z,- 


•000 C3^,- 


•14c2^^ = 


±.023, 


(124d) 


■'^^ '^IZZ ~ 


- .31 C2ZZ ± 


.52 c^zz ~^ ■ 


OAQC2Z,- 


•000C3^,± 


•35 C2^^ = 


±.079, 


(124e) 


.023 c,zz^ 


- .88 C2^^ ± 


.19 c-^zz - ■ 


092 C2Z,- 


•000C3^,± 


•43 C2^^ = 


±.095. 


(124f) 
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All the errors in eq. (11240 are approximately a factor of 2 larger than the corresponding ones in 
eg. (11221) . mainly because of the loss of information due to the binning and also to the reduction 
of the dimensionality of the weight functions from 4 (xj, cos6'i, cos 6*2, 0) to 2 |t2|)- 

3. Sensitivity at ^/s = 500 GeV and 1 TeV 



At = 500 GeV with L = 100 fb"^ with \P\ = 80%, we find : 



002 + .000 C2^^ - .007 C3^^ - .004 C2^^ + . 


99995 c^z^ - 


-.005c2^^ = 


± .00057, 


(125a) 


.034 q^^ 


- -20 C2ZZ - -064 c^zz - -61 c^z^ 


-.mic.z. 


- .76 C2^^ = 


±.0029, 


(125b) 


.044 


- -28 c^zz - -030 c^zz + -78 c^z^ 


+ .000c3^^ 


— .55 q^y-y = 


±.0035, 


(125c) 


.19 c^zz ^ 


- -31 C2ZZ - -93 c^zz + -045 c^z^ - 


-.007c3^,- 


-■029c2^.y = 


±.0067, 


(125d) 


.34c^zz 


- .84 C2ZZ - -23 c^zz - -098 c^z^ 


-.000c3^, 


+ .33 q^y-y = 


±.013, 


(125e) 


.92 c^zz^ 


- .27 C2ZZ + -28 c^zz + -012 C2z^ - 


-.000c3^^- 


- .063 C2^.y = 


±.031. 


(125f) 



The cross section now is 6.6 fb, or 660 events with 100 fb~^; see figEl We therefore expect that 
the statistical errors should be about 2/3 of those at ^/s = 350 GeV. In fact, all the errors of 
eq. (11251) are smaller than the corresponding errors in eq. (11221) by a factor of 1.2 to 2.4. At this 
energy, as mentioned above, kinematic difference in the Z and 7 propagators diminishes, and 
the H'-f'-f coupling 03^^ tends to mix with C2z^ and the HZZ couplings. 
And at y/s = l TeV for L = 100 fb"^ with \P\ = 80%, we find : 



002 c,zz + -002 C2ZZ - -020 c,zz " -019 c^z-y + -9996 c,z^ - 


-.006c2^^ = 


± 


.0004, 


(126a) 


.008 c^zz + -033 C2ZZ - -15 c^zz - -95 c^z^ - .022 c^z^ 


-.26c2^^ = 


± 


.0018, 


(126b) 


.089 c^zz + ■^'^c^zz - -^Tc^zz + ■'^^c^z^ - -017 03^^ - 


- .09102^^ = 


± 


.0022, 


(126c) 


.035 - -Q^c^zz + ■O'^^c^zz + -1702^^ ± -000 Cg^^ 


-•72C2,, = 


± 


.0025, 


(126d) 


.15 Cizz ~ -72 C2^^ — .17c3^2 — .17 C2^^ — .002 C32^ 


±-64c2^^ = 


± 


.0058, 


(126e) 


.98 c,zz + -12 C2ZZ ± -11 c;zz ± -OH c^z^ + -000 c^z^ - 


- .063 C2^^ = 


± 


.019. 


(126f) 



The cross section is around 13 fb, or 1300 events with 100 fb~^. The error of c^Z'y eq. (ll26al) 
decreases by a factor 1.4 from the corresponding one in eq. (ll25al) for 500 GeV. The other errors 
are reduced more significantly. At high energies, the distinction among the weight functions 
becomes clear and c^z-y^ ^2^75 '^szz ^izz measured rather independently, while only C2zz 
and Cg^^ are correlated. Since the 7 and Z propagator factors are similar at high energies, the 
distinction come from the extra (^1,^2) dependence in eq. (lll7p and from the azimuthal angle 
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dependence in the factor Q and Q'; see eq.f lBSP in appendix B. Most significant improvement 
is found for the error of the eigenvector with dominant c^zz component, which decreases by a 
factor of 3 from 0.0067 at = 500 GeV in eq. (fT25dl) to 0.0022 at = 1 TeV in eq. (fT2fe|l . 

Overall, we find that the HZ'y couplings are most accurately measured at high energies 
in the double-tag eeH process. However, the accuracies of the Cg^^ and Cg^^ measurements 
cannot compete with those from the ZH process, as can be seen e.g. by comparing eqs. fll26al) 
and (]126bl) with eqs. fllOQal) and fllOQcl) . respectively, at = 1 TeV. The double-tag events 
are found to be most important in distinguishing c-^^z^ ^^^^ HZZ coupling, from C2zz 

and c^zz- 



VII. NO-TAG ee ^ {ee)H PROCESS 

Although the double-tag t-channel eeH process shows certain sensitivity to the H'-f-f coupling 
Cg^^, the sensitivity at y/s = 500 GeV and 1 TeV with Lq = 100 fb~^ is not sufficient for resolving 
the if 77 coupling of the SM from the W^-boson and top quark loops. We can expect better 
sensitivity to the H'j'j coupling when the are not tagged, because of the fusion of almost 
real photons. The produced Higgs-boson has little px, and can be distinguished from the UeU^H 
events. 



A. Cross section 

The momentum and helicity assignment for the initial and final particles in no-tag {ee)H 
process is the same as the double-tag eeH process, eq. flllOp . The momentum and helicity of 
the intermediate almost real photon is assigned as 







Ai 


-) 

2 / 




2 






As 


2 ) 




2 



ki,^) ^ e-(pi,^)+7(fc'i,A;), (127a) 

-l{k',,K). (127b) 



where A'^ and Ag is the helicity of the equivalent real photon emitted from e and e"*", respectively. 
The final electron and positron escape detection by going into the beam pipe, and an almost real 



photon of the virtuality <^ rn^j is emitted in the same direction; see e.g. ref. 26|] for a review. 
The Z boson contribution is suppressed by a small factor of |g^|/m| <^ 1. The minimum and 
maximum magnitude of the square of the electron or positron momentum transfer is calculated 
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as 



ymin 



^imax 



z) - [1 — cos6'r, 



'128) 



l-z' ~ ' 2 

where z = 2 E^j ^/s is the energy fraction of the photon, and 6'min is the polar angle below which 
the final escapes detection into the beam pipe. In our analysis, we set cos 6 
eq. fl33|) . The equivalent real photon distribution can be split into two parts 



D, 



7/e=l 



^^'^^^'z,(f) + ir^ 



where the helicity preserving component {aX' 



are 



27| 



,/e± v^,Q'), (129) 
and the helicity flip component {aX' = — ) 
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2T 

a 



In 



Zf 



In 



l]+z 



- 1 



^min 



respectively, with 



Q2 = min{m|^, q^^^J. 



(130a) 
(130b) 

(131) 



It is worth noting that the last term z in the parenthesis of the helicity preserving (crA' = +) 
distribution f ll30b ) comes from the helicity flip {aX = — ) amplitude where the final helicity 
is opposite from the initial one. Since contribution from the the momentum transfer much 
smaller than 1 GeV^ dominates the photon distribution, we use a = 1/137 in eq. fll30p . The 
real photon approximation for the matrix elements cease to be valid when = \ti\ becomes 
0{'m'jj). In our study, the maximum momentum transfer g^^^ stays below ~ 0.25 m'jj even at 
^/s = 1 TeV, and hence our predictions do not change significantly when we replace the scale 
from by 0.1 m^, in eq. fll3ip . 

The cross section for the no-tag {ee)H events is now expressed as 



a 



{ee)H 



dxi [ dx2 D'^2~{xi,Q'') D 
Jo 



where 



16 TT' 



77 



m 



6{s 



m 



H) 



(132) 



(133) 



H 



Here, Xi and are the energy fractions and the helicities of the equivalent real photon in 
e~{i = 1) and e"^(i = 2), respectively, a/2 and (3/2 are the e~ and helicities, and s = X1X2S 
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is the square of the 77 colhsion energy. Higgs-boson is produced only when X[ = A2. Therefore, 
the total cross section and the rapidity i/h = 1/2 log(a^i/a^2) distribution of the Higgs-boson, 
which has strong dependence on the helicity profile of the equivalent real photon distribution, 
depends on the beam polarizations P and P. We find, however, that the polarization effects 
survive only when both e~ and have non-zero polarization. 

The un-polarized cross section of the no-tag {ee)H events in figJS] for = 120 GeV and 
those of figJUfor rrijj = 100 GeV to 200 GeV are calculated by using the above formula with 

r,, = l^r. (134) 

256 71-^ ' ' 

where the loop function / sums over the contributions from the VT-boson and the top-quark 



291]. The enhancement due to the W^-boson pair production threshold can be observed in figJH 



For = 120 GeV, the decay width of eq. fll34p gives = 0.0073 MeV, and the cross section 
is 0.15 fb, 0.32 fb, 0.56 fb and 1.3 fb at = 250, 350, 500 and 1000 GeV, respectively, as 
shown in figlSl When the integrated luminosity is 100 fb~^, the expected number of events are 
15, 32, 56 and 128, respectively. Therefore, in this study, we consider the no-tag events only at 
^ = 500 GeV and 1 TeV. 



B. Error Estimation for C2^^ 

The interaction Lagrangian for the effective if77 vertex 

^ii2^ = ~ "^277 H A^y A^^ , (135) 
contains both the SM loop contribution and new physics effects. We parametrize them as 

c^,, = cff, + /^c^,,. (136) 

The SM contribution is approximated by the loop function /, and its numerical value for = 
120 GeV is 

C2^'7 = ^ ^ = - 0.00187, (137) 

with negligibly small imaginary part from the b and r contributions. Our approximation of 
the point-like coupling in the effective Lagrangian eq. (11351) is valid as long as ^ m^,mj. 
The cross section of the no-tag {ee)H process for rrifj = 120 GeV is expressed in terms of the 



49 



generalized if77 coupling given in eqs. (11351) and fll36p as 

^ee^(ee)H^y^ = 500 GeV) = 1.61 X 10^ (^4^^ + Ac^^^^ 

^ee^iee)H^^ = 1 TeV) = 3.70 X 10^ (^4^^ + Ac^^^^ 



fb. 



fb. 



(138a) 
(138b) 



As mentioned earlier, this cross section does not depend on the e~ polarization P, as long as 
the polarization is zero, P = 0. 

While evaluating the error of the cross section, we take into account of the contribution 
from H^iy-fusion Ve^eH events, which cannot be distinguished from the no-tag {ee)H events 
when the Higgs boson has small transverse momentum. In the following analysis, we make a 
very naive estimate for the background and the errors by assuming that the Higgs boson 
can be resolved at 3 GeV accuracy and that the SM background contribution can be estimated 
from the high distribution. In figJHl we show the expected p^ distribution of the Higgs 
boson from the no-tag {ee)H process as the red-solid rectangle, which is overlaid above the 
H^VT-fusion contribution shown by the gray-shaded histograms; figlH^ is for ^/s = 500 GeV and 
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FIG. 8: The histograms for the number of events at low Higgs px 



figHb is for y/s = 1 TeV. Also shown in the figures are the statistical errors corresponding to 
the background lyVF- fusion events, and that of the sum of the no-tag {ee)H and lypl^-fusion 
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events for the lowest pt (pt < 3 GeV) bin with Lq = 500 fb ^. We can then estimate the error 
of the signal {ee)H cross section as follows : 



0"sig — CTtot — O'BG — O'tot — O'BG ± 











iVtot 


iVBG 



= asig±y , (139) 

where a represents the mean value and L is the integrated luminosity. The errors of the total 
and the background cross sections in the lowest bin are estimated to be 0.023 fb and 0.012 
fb, respectively, at y/s = 500 GeV, and 0.033 fb and 0.016 fb, respectively, at i/i = 1 TeV. By 
taking the squared sum of the errors, we obtain the following estimates : 



(J: 



sig 



cr.. 



Slg 



0.56 ± 0.045 fb at = 500 GeV; L = 500 fb~\ (140a) 
1.29 ± 0.064 fb at = 1 TeV; L = 500 fb"\ (140b) 



We then find from eg. (11381) 

cfi + Aco^^ = - 0.00187 j ^^-^^^^^^ at = 500 GeV; L = 500 fb"\ (141a) 

I -0.000077 

cfi + Ac. = - 0.00187 j +'^-°'^°°^^ at y/s = l TeV; L = 500 fb"\ (141b) 

2t7 277 I -0.000047 V , y ; 

The above results show that the no-tag {ee)H process at e^e~ collider can measure the H'^'y 
coupling with the 4% and 2% accuracy with the integrated luminosity of 500 fb~^ at y/s = 
500 GeV and 1 TeV, respectively. In contrast, the sensitivity of the double tag eeH events on 
the C2y-y coupling in eqs. (11251) and (11261) . gives the error of 0.0016 and 0.0014 at y/s = 500 GeV 
and 1 TeV, respectively, with 500 fb~^, which is barely enough to resolve the SM coupling at 
1 a level. From this exercise, we can conclude that the no-tag {ee)H process is about 20 to 30 
times more sensitive to the H'-ff coupling than double-tag eeH process, at ^/s = 500 GeV and 
1 TeV, if our naive estimation of the background is valid. 

This naive estimation of the errors is based on the assumption that the background con- 
tribution from the VTVT-fusion process can be estimated from the high distribution. This 
assumption can be tested experimentally by using the different polarization dependence of the 
signal and the background. The background W^VT-fusion process scales as 1 — P, i.e. , it in- 
creases for P < and decreases for P > 0, while the signal remains independent on P as long as 
positron is unpolarized. Once the positron polarization is available, both the total cross section 
and the rapidity distribution of the no-tag {ee)H events depend on the sign of the product 
P ■ P, which should be a distinctive signature of the 77 fusion. 
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Although the sensitivity to the if 77 couphng is expected to be measured far more accurately 
with the photon-linear collider 30|, ISl], our results show the sensitivity limit when it is not 
realized. 



VIII. SINGLE-TAG e+e" {e^)e^H PROCESS 

This process essentially measures e^7 — e^H , where the initial photon comes from the 
other e^ beam as an equivalent real photon. Since we have studied both double-tag eeH and 
no-tag {ee)H events, it is worthwhile to examine the single-tag {e)eH process for completeness. 




e [ki) or e+(fc2) e {pi) or e+(p2) 

FIG. 9: Feynman diagram of 6^7 — > e^i? 



A. Cross section 



The momentum and the helicity assignments of the single-tag {e)eH process is also the 
same as those in the double-tag eeH process in eq. flllOl) . We use the equivalent real photon 
approximation for the current from the untagged e^. Since the helicity of the equivalent real 
photon depends on the parent helicity, and also the electroweak interactions of the tagged 

current are sensitive to their polarization, we give the general expression of the cross section 
with both e~ and e+ beam polarizations, P and P, respectively; 



(e)eH 
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a,f3 



l + aP\ fl + f3P 
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e 7- 



7/e 



6^7- 



7/e 



(142) 
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Here the polarization-dependent real photon distributions are as defined in eq. fll30p . and the 



e^7 e^H cross sections are 



1 

2^ 



OA' I ^ 



32 TT s 



1 - 



m 



H 



d cos 9* 



(143) 



where s = s a/ 2 is the helicity, A'j is the photon helicity, and 6'*± is the polar angle of 
the final tagged in the colliding e^7 cm. frame, measured from the incoming momentum 
direction. The integration region of cos is constrained by the tagging criteria of eq. (132!) 



COS^pi < 0.995 = COS^ri 



(144) 



in the e^e collision cm. frame. We find 

~ 1 + T r. r-r. n — ttttt < cos 9„± < 1 

[2;+ (1 -2) (1 +C0S^„,in)/2] 

The interaction Lagrangian relevant for the e^7 

1 . .„., 2 , 



1 — cos 6„ 



z+ 1- z) 



e^H process is 



cos 6'min)/ 2] 



• (145) 



V V 



(146) 



The coupling C2^^ dictates the 7-exchange amplitude and the combination 03^^ — 03^^ dictates 
the Z-exchange amplitude in figlH The helicity amplitudes for the e~7 — > e~H process are 
expressed as : 
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\^2Z-y ^3Z-f) 



sm ■ 



H 



COS 9* 



is the massless t-channel current of 
eq. flll3ap . ti = {ki — piY = —s Xi (1 — cos6'*)/ 2, xi is the energy fraction of the tagged e~, 
9* is the scattering angle in the e~7 rest frame, ^2 is the four momentum of the equivalent real 
photon emitted from e"*". The amplitudes for e"''7 — e~^H are obtained from eq. fll47p simply by 



(147) 



where the T^^^ form factors are as in eq. ffT^ . J^-^-^^-) 



replacing a, 



Zee 



ti = {h - Pi? and A'2 by /3, h 



(^2 ~ P2Y and A'^, respectively, where 



9* is now the scattering angle between the initial and final e+ in the €"''7 rest frame. 

It is instructive to show that the single-tag {e)eH cross section of eq. fll42p reduces to the 
no-tag {ee)H cross section of eq. fll32p in the limit of ^ m|. In this limit, the 7-exchange 
contribution dominates the amplitude, and by noting cos 6** ~ 1 and xi = 1 — rrijj/ s, eq. (11431) 
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becomes 
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(148) 



Comparing the final expression of eg. (11481) witli eg. (1130b ). we find that the term proportional 
to in the uA = + distribution from the electron helicity flip amplitudes, is not reproduced, 
since we have neglected the electron mass of the tagged e~. Inserting eg. (ll48p in eg. (11431) . we 
reproduce the cross section (11321) for the no-tag {ee)H events. 



B. Error estimation for C2^^ and Cg^^ — c^z^y 

To calculate the SM 1-loop contribution, we again adopt the approximation of replacing the 
1-loop vertices by localized effective couphngs, 

%7 = <7 + ^%7' ^2z, = + Ac,zr (149) 

It is important to note that since the virtuahty of the tagged currents (t) can be larger than 
the weak boson mass scale that dictates the spatial extension of the vertex, the approximation 
of eg. (11491) is not as excellent as that in eg. (11361) for the no-tag {ee)H events. In evaluating the 



SM contribution, we therefore retain the t-dependence of the loop function 321], 



and use the t-dependent 'couplings' in the evaluation of the 67 eH matrix elements. In 
other words, only the new physics contributions via Ac2^^ and Ac2^^ are assumed to be local 
in eg. (ll49p . With regard to the if 77 couplings, C2^(0) = —0.00187 in eg. (11381) determines the 
if ^ 77 width, while it reduces to C2^(— m|) = —0.00160. Similarly, for the HZ'-f coupling, 
clf^{ml) = -0.00322 determines the H ^ Z-f width, while it reduces to c|f^(0) = -0.00276 
and c^^{—my) = —0.00243, all for = 120 GeV. As a conseguence of the W boson 
dominance in the loop, we find that the ratio 

C2Z7(^I)/ c|^(0) - 1-7, is essentially the 
ratio of the W boson gauge couplings gzww/g-yww — 1-8- Since the effective couplings are 
gauge-dependent when they are away from their on-shell limit, t = for CgM, and t = m| for 



C2z\, respectively, the SM cross section based on the t-dependent effective couplings should 
be regarded only as an order of magnitude estimate. However, the box diagram contribution 
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to the 67 eH process, which is necessary to recover the gauge invariance, is found to be 



neghgibly small numerically 32 1. 



Because the single-tag {e)eH process has small cross sections as shown in figJS] and figJH 
we use only the total cross section to constrain the H'-ff and HZ'y couplings. In order to 
quantify the effects of the e~ beam polarization, we obtain the following parametrization for 
the single-tag e~ events for the e~ helicity a and un-polarized e+'A; 



Here the average over e"*" helicities is replaced by the sum over the photon helicity A'^. We find 
for = 500 GeV, 

an = 15000 {c^.^Y + 890 {c^z^ - c^^.f " 5500 c^,, {c^z, - c^z,) = 0-025 fb, (152a) 
<JL = 15000 {c^^^f + 1200 {c^z^ - c.,z^f + 6400 c^^^ {c^z^ - c^^^) = 0.066 fb, (152b) 

and for = 1 TeV, 

(jR = 24000 {c^^^f + 2900 {c^z^ - c^^^f - 15000 c^^^ {c^z^ - c^^^) = 0.022 fb, (153a) 
aL = 24000 {C2,,f + 3900 {c^z, - Hz^f + 17000 c^^^ {c^z, - c^z,) = 0-107 fb. (153b) 

The large coefficients of (cg^^)^ indicate that the 7-exchange amplitude dominates the cross 
section even for the tagged events. The effects of the HZ'-y couplings can be measured from 
the interference between the 7-exchange and the Z-exchange amplitudes, which is destructive 
for cr and constructive for e^, reflecting the relative sign of their couplings, 9]^^= q'r^ = — e 
and f^f^"^ = (—1/2 + sin^6'w), gft^ = 9z sin^6'w- Therefore, we expect that the beam 
polarization can be a powerful tool to distinguish the HZ^ coupling from the iir77 coupling. 

It is also worth noting that the cross section after the e tagging condition of eg. (11451) can 
become smaller at high energies, as we find for gr in eq. fll52b ) and eg. (1153b ). In figJTOl we 
show the rapidity distribution of tagged e~ in the laboratory frame. Although the total cross 
section increases monotonically with -y/s, the distribution shifts to large r/ as ^Js grows. It 



[t] The coefficient Aa is determined by evaluating the 7-cxchangc contribution only, such that the SM cross 
section is reproduced with a constant Cj^^, which may be interpreted as the average of (c2^^)^ = (c2^(t)^)77- 
Likewise, the coefficient Cq is determined by requiring that the SM contribution to the 7-Z exchange interfer- 
ence is reproduced for a constant value of c^^^ (^2^7 " ^zz-^ ={'^^■^{^^1 iik'^^-^^)) z-i ■ Finally, the coefficient 
Ba is fixed by demanding that the total SM cross section is reproduced for the above average values. 
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FIG. 10: The rapidity distributions of tagged e in the lab. frame for each initial helicity state and at 
each collision energy. The thick solid line shows the detection limit, above which e~ escape detection. 



implies that more and more fraction of events escape detection (cos^e > 0.995 or rje > 2.99) at 
high energies. The suppression of the cross section for (?7e < 2.99) as shown in fig{TUl can be 
attributed to the cancellation between the 7- and Z- amplitudes for the tagged e^, with large 
scattering angle in the e~7 rest frame. 

In terms of the above parameterizations eqs. fll52p and fll53p . we can determine the total 
cross section for polarized e~ beam and unpolarized e"*" beam as follows : 

P = 0) = ^ (an + ctl] . (154b) 



2 

At = 500 GeV with L = 500 fb~\ we find 



^e-H^p ^ p ^ 0) = 0.030 (± 0.011) fb, (155a) 

^e-H(^p ^ _Q_g^ P = 0) = 0.062 (± 0.016) fb, (155b) 
^e+H(^p ^ P = 0) = 0.046 (± 0.0096) fb, (155c) 



and at y/s = 1 TeV with L = 500 fb ^, we find 



a"' 


"^(P = 


+ 0.8 


,P 


= 0) 


a'' 


"^(P = 


-0.8 


,P 


= 0) 




"^(P = 


±0.8 


,P 


= 0) 



0.030 (±0.011) fb, (156a) 
0.098 (±0.020) fb, (156b) 
0.064 (±0.011) fb. (156c) 



56 



The errors in the parentheses are for L = 250 fb ^ each for P = 0.8 and P = — 0.8. 
From eqs. fll52p and (11551) . we find 



C2Z7 - C3Z7 = - 0.00201 ± 0.0012 / 1. 
C2^^ = -0.00166 ± 0.00016 \ - .71 1. 



(157) 



at ^/s = 500 GeV. The mean values may be interpreted as the average of the t-dependent SM 
couphngs in the tagged events, which turn out to be approximately c||'^(— m|) and c|^(— m|), 
but with shghtly smaller magnitudes. The smallness of the magnitude of (c|^^(t)) reflects the 
larger mean value of the momentum transfer in the 7-Z interference than in the purely 
7-exchange contribution. When we compare the result eq. (ll57p with that of the no-tag {ee)H 
events in eq. (ll41al) . we find that the no-tag events are about a factor 2 more sensitive to the 
H'-f'-f coupling. On the other hand, the sensitivity to the HZ'-y coupling, ±0.0012, is one order 
of magnitude worse than that of the ZH process in eq. (fTUHj) rescaled for L = 500 fb~\ Hence 
the measurements of the single-tag {e)eH events do not improve the constraint on the HZj 
coupling significantly. 

For -y/i = 1 TeV, we find 



- C3Z7 = - 0.00180 ± 0.00063 
C2^^ = - 0.00145 ± 0.00014 




'158) 



from eqs. (11531) and (11561) . The magnitudes of the mean values decrease slightly from those of 
eq. (11571) at ^/s = 500 GeV, because of the higher value of the typical momentum transfer, \t\, 
at ^/s = 1 TeV. The error of Cg^^ is almost a factor of 3 larger than that of the no-tag {ee)H 
events, eq. (1141bl) . while that of Cg^^ is more than one order of magnitude larger than that from 
ZH events in eq.( 1109p . rescaled for L = 500 fb~^. 

Although the process 67 — > eH gives us clean measurements of the HZ'-y couplings with 
e~ beam polarization, the sensitivity is rather low in e^e~ collisions because of the small and 
soft photon flux from the bremsstrahlung. This process will become more important once the 
photon linear collider option is realized 32 1. 



IX. LUMINOSITY UNCERTAINTY 



So far, the errors and their correlations based on the optimal observables method for Higgs- 
gauge boson effective couplings are computed by assuming the true luminosity L. However, the 
error in the measurement of the L can affect the measurements of some effective couplings. We 
attempt here to study the impact of the luminosity uncertainty on the precision measurements 
of the HVV couplings. 
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In the presence of the luminosity uncertainty, the true luminosity L can be estimated as 

L = fL, / = 1±A/, (159) 

where L is the measured mean value, and A/ is its 1-a uncertainty. The function given in 
eq. (|2^ of section UTTl should then be redefined as follows : 
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is now {n + 1) X [n + 1) matrix with 

Co = / - 1 



Ssm($) 



So(<f ) = SsM(<f ). 



, (160b) 

(160c) 
(160d) 

(161) 
(162) 



It is straightforward to integrate out the Cq = / — 1 dependence and obtain the probability 
distribution of the parameters ci to c„ in the presence of the luminosity uncertainty. In our 
study, we note that the two couplings, c^yyi^ and c^^^ have the weight functions that are 
identical to the SM distribution. Because of this, we can study the impacts of the luminosity 
uncertainty algebraically by using the functions written in terms of c^^^^ and c-^^zz- 



^^...w-(*) = 2Ssm($), 
for the lyW^-fusion process given in eg. (1451) . and 

S,^^^(<l>) = 2SgM($), 



(163) 



(164) 



[I] For no-tag {ee)H and single-tag {e)eH processes, the H^^ and HZ^ couplings of the SM at one loop order 
are used to calculate the SM cross-sections. There, however, the statistical error is dominated, and the errors 
due to the luminosity uncertainty can be safely neglected. 
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for the ZH production and double-tag eeH processes given in eq. (173|) and eq. flllOp . respectively. 
We can express the total function with the luminosity uncertainty given in eg. (11611) for a 
particular cm. energy, i/i as 

HWW 
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2/ ' _ 
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\ ^"^ 1 




L(A/)v/iJ 



(165) 



In the following analysis, we assume that the luminosity uncertainty is common to all the 
processes at each collision energy -y/s. 

Let us examine the effects at ^/s = 500 GeV in some detail. The function for the HWW 
process in eq. (l62|) can be expressed in the form 
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(166) 



where c'^vviy = '^iww + (/ ^ 1)/ 2, c'^Y/^^ = Ciww {i !)• Likewise, the function from the 
ZH and the double-tag eeH processes in sections |V] and |Vl] can be expressed as 
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2 / / 

X I Cizz ~^ (^izz 
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(167) 



kl 



where c[z z = Ciz z + if -'^)/ 2 and 4, 4, 4, ci,, 4= c^zz^ c^zz^ c^zv ^327. %7- Now, the 
luminosity uncertainty in the x^ function of eg. (11651) at ^/s = 500 GeV can easily be factored 
out as 
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(168) 

(169a) 
(169b) 



k=l 



and the reduced function is 

x'(q)500 = x'(q)500 - [(A/)^50o]'^'- 



(170) 



We can use the reduced function to study the constraints on the non-standard couplings in 
the presence of the luminosity uncertainty. It should be noted that because of the last term 
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of eq. fll70p there appear correlations between the HWW couphngs and the HZZ/HZ-f/H'yj 
couphng measurements. Therefore, the effects of the luminosity uncertainty is important when 
we study the constraints on the coefficients of the dimension-6 operators, eq. ffT^ . since some 
of them contribute to both couplings. 

The reduced function of eg. (11701) gives a 9 x 9 covariance matrix. The results for ^/s = 
500 GeV, L = 500 fb^^ and A/ = 0.01 can be expressed as, 

/ 1 
.990 1 

.991 .995 1 

.043 0. 0. 1 

0. 0. -.63 1 

0. 0. -.81 .033 

0. 0. -.006 -.063 

0. 0. .0 .050 

0. 0. .22 -.27 



.0 
.0 
.0 
.0 
.0 
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.043 

-.046 

-.28 



-.83 1 
-.040 .034 1 



(171) 



When the above results are compared with those without the luminosity uncertainty in eg. (1631) 
for the HWW couplings and the corresponding expression for the HZZ/HZ'-f/H'-f'-f couplings, 
after rescaling of the errors by l/^/b ^ 0.447, we find the followings. The errors of c^t^^^ 
and c^zz slightly larger than 1 / \/h of the corresponding statistical errors, while all the 
other errors are not affected much. The correlations between c^^-^ and Cg^j/iy or Cg^^^y are 
reduced slightly in magnitude from eg. (1631) . and those correlations between c^zz the other 
HZZ, HZ^, H'-f'-f couplings are also reduced slightly in magnitude. Finally, the off-diagonal 
sub-correlation matrix in eg. (ll71|) is almost vacant except for the {cnyw^'^izz) component, 
which shows positive correlation because the luminosity uncertainty affects both Cmr^/ and 
c^zz couplings in the same way. All the effects of the luminosity uncertainty are rather small 
because the statistical errors of c^^y^^ and c^zz ±0.045 and ±0.013, respectively, which are 
significantly larger than the error from the postulated luminosity uncertainty, A//2 = 0.005. 
When the errors due to the luminosity uncertainty become dominant, the errors of Ciyy^^^ and 
c^zz '^ill stop decreasing with the luminosity, and the correlation between c^^y^ and c^zz ^^^^ 
grow. 

It is instructive to study the impact of the luminosity uncertainty analytically in a very 
simplified example where only the two couplings, c^^yy and c-^zz^ retained in the amplitudes 
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at ^/s = 500 GeV. In this limit, the reduced function of eq. fll70p is simply 

X^(Ci)500 = clwW [(^SCkT^) ]i1+CizZ [{^500 ^) ]ii 
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where the combined error of (A/)^||p is 
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(172) 
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As noted above, the second term in eg. (11 72 p generates the correlation between the errors of 
^iww ^^"^ ^izz- -^o^ example, if we take the limit where the statistical errors are much smaller 
than the luminosity uncertainty, {Acmry/Y, {Ac^zz?' {^fY i ^^e above reduced function 
can be expressed as 
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• (174) 



In the leading term, only the combination {c\\/i/^ji/ — C\zz) constrained, and the next-to-leading 
term proportional to 1/ (A/)^ constrain the combination (Ac^^^^)^ (^iww + (Acivi/h^)^ ^^izz- 
In the limit of large statistics, the first term dominates the function and the correlation 
approaches the unity. 

In figHH we show the errors of c^^^^ and c^zz couplings as a function of the integrated 
luminosity L for the luminosity uncertainty A/ = 0.01 at = 500 GeV. The solid line shows 
the error of the combination Cy^^ — c^zz which decreases as 1/ asymptotically, since the 
luminosity uncertainty is canceled out by taking the difference. On the other hand, the dashed 
line for Cj^jy^j/, and the dash-dot line for c^zzi approaches the same value (A/)/ 2, showing 
that the individual errors will be dominated by the luminosity uncertainty. More explicitly, the 
errors of the couplings, c^y^^ and c^zzi '^^^ ^e expressed as 



^YWW — =1= 
with the correlation 
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FIG. 11: Error of c^y^y^ and c^zz integrated Luminosity at ^ = 500 GeV when the luminosity 
uncertainty is A/ = 0.01. 



X. CONCLUSION AND DISCUSSION 



The imprint of the dynamics of the symmetry breaking physics is inherent in the interactions 
of the Higgs boson and the gauge bosons. In this article, we attempt to evaluate the potential 
of the future linear e^e~ colliders, such as the ILC, in probing the dynamics of all the CP-even 
and gauge-invariant dimension six operators of the Standard Model (SM) fields that affect the 
Higgs-gauge-boson couplings when there is one SM-like light Higgs-boson. For this purpose, 
we study all the processes that are sensitive to the HW couplings: the z/gZ/g-ff production via 
lyiy-fusion is sensitive to HWW couplings, the ZH production process is sensitive to HZZ 
and HZ'-f couplings, the double-tag eeH production process via t-channel Z and 7 exchange is 
also sensitive to HZZ, HZ'y and H'-ff coupling, the no-tag {ee)H process measures the H'yy 
coupling, and the single-tag {e)eH process measures HZ'-f and if 77 couplings. In order to 
quantify the resolving power of each process, we allow all the effective HVV couplings to vary 
freely in the fit and adopt the optimal observables method to constrain them in each process, 
at a few selected collision energies (a/s =250, 350, 500, 1000 GeV), and with or without e~ 
beam polarization. 
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A. Summary of the constraints on the HVV coupUngs 



Here we summarize our results for the effective HVV couplings at the ILC. 

• All our results have been presented for nominal integrated luminosity of Lq = 100 fb~^, 
except for the no-tag {ee)H and single-tag {e)eH events, where we give our error estimates 
for L = 500 fb~^ at = 500 GeV and 1 TeV only. 

• Most of our results have been given for e~ beam polarization of |P| = 80% with no 
e"*" beam polarization |P|=0, where exactly half of the total integrated luminosity L is 
delivered with P = \P\ and P = —\P\. 



42], but presented 



• Results for \P\ = \P\ = have also been calculated for all the cases 
only for the ZH production process at y/s = 250 GeV; see eq. fll04p . 

• Our results for the HWW effective couplings are shown in eq.( l66i) for y/s = 250 GeV, 
eq.dSZD for = 350 GeV, eq.([62]) for = 500 GeV and eq.([68D for = 1 TeV. 

• The results for the HZZ and HZ'-/ couplings in ZH production process with |-P| = 80% 
are shown in eq. ffT06|) for = 250 GeV, eq.ffTOTD for = 350 GeV, eq. ffT08|) for 
^ = 500 GeV, and eq.(^M) for = 1 TeV. 

• The results for the HZZ, HZ'-/ and H'-f-f couplings in double-tag eeH process from t- 
channel vector boson fusion with \P\ = 80% are shown in eq.( ll23p for ^/s = 250 GeV, 
eq. ffTMj) for ^ = 350 GeV, eq.(^2B for = 500 GeV, and eq.(^2B for = 1 TeV. 

• The H'j'y coupling results from the no-tag {ee)H process are given for L = 500 fb~^ in 
eq. ffrm at = 500 GeV, and in eq. ffTilbll at = 1 TeV. 

• The HZ J and H'j'^ coupling results from the single-tag {e)eH process for L = 500 fb~^ 
with \P\ = 80% are given in eq. ffTHTD at y/s = 500 GeV, and in eq. ffT58ll at ^/s = 1 TeV. 

• The impact of the luminosity uncertainty, L = fL with / = 1 ± A/, should be taken into 
account according to eq. (ll65p . where the coefficients Cmriy and Ci^z the function at 
each energy is replaced by c'^yy = c^yy+{f—l)/Af and by adding the term (/— 1)^/ (A/)^ 
at each energy. After squaring out the (/ — 1) dependence, we obtain the reduced 
function at each energy. 

The combined analysis in this section are performed by gathering all the above results, by 
adding all the functions with appropriate weights. 
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In particular, we report the following two cases; 



2/ N -^250 ~2, \ , -^350 ~2, \ 

X [Ci)i = — — X (Ci)250 + X (Ci)350 

2( \ 2( \ , -^1000 

X (q)i+ii = X {ci)i + — — 



^0 

X^(ci)iooO) 



-^500 ~2 / N 

^— X (Cijsoo 



(176a) 
(176b) 



where x^(ci)i gives the combined results of the ILC phase I with the maximum energy of 
= 500 GeV and x^(cj)i+ii gives that of combining the results from ILC-I and ILC-II at 
^/s = 1 TeV. The reduced function at each energy, x^(cj)^, are obtained for the luminosity 
uncertainty of (A/)^ at each energy. 

Since we present the individual covariance matrix separately in each process at y/s = 
250 GeV, 350 GeV, 500 GeV and 1 TeV, in terms of its eigenvectors and square root of 
eigenvalues, we can estimate the constraints for an arbitrary integrated luminosity with an 
arbitrary luminosity uncertainty A/ at each energy. As an example, we show the results for 
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when the luminosity uncertainty is A/ = 1% at all energies. 

Using the function defined in eq.f ll76ai) . we evaluate the errors and correlations of the 
effective couplings for ILC-I, with the integrated luminosities given in eq. (11770 and obtain 
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When we compare the above combined results from the three ILC-I energies with those of 
eq. (11711) at ^/s = 500 GeV, we observe the foUowings. As for the HWW couplings, we observe 
reduction of the errors of all the 3 couplings and that of correlations. This is essentially because 
contributions have different dependence on ^/s. As for the HZZ/HZ'-f/H'-f'-f couplings, the 
errors of Cg^^, c^zz^ ^2Z-y C3Z7 essentially determined by the ^/s = 500 GeV experiment, 
while that of c^^^ is reduced significantly by including the lower energy data, because the ZH 
production cross section is larger at lower energies, see figOl The error of Cg^^ is essentially 
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determined by the no-tag {ee)H events.We find that the single-tag {e)eH events, eq. fll57p . do 
contribute significantly to improve the HZ^ couplings measurement. 

Only the errors of c^w^y and c^^^ and the corresponding rows of the correlation matrix 
are affected by the luminosity uncertainty A/ = 0.01 at each energies, assumed in the fit. 
By combining the results of eqs. fl64l) . (I^Sj) and fl62l) we obtain the combined error of c^-^y^ 
without the luminosity uncertainty to be 0.025, which is still 5 times larger than the error due 
to the luminosity uncertainty, A//2 = 0.005 at each energy. Therefore, the error of Cj^^^/^y 
does not increase significantly with the inclusion of the luminosity error. On the other hand, 
by combining the results at ^/s = 250, 350 and 500 GeV, we find that the combined error of 
'^izz without the luminosity uncertainty is 0.0045, which is comparable to the error due to the 
luminosity uncertainty. We therefore find almost 30% larger error for c^^^ eq. fll78p . The 
correlation between c^^r^r and c^zz ^'^^ 0.17, which is still small because the statistical error 
of 

Ciww is much larger than the error due to the luminosity uncertainty. 
After combining the ILC-I and ILC-II results with the integrated luminosities of eq. fll77p . 
we find 
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The HWW couplings are measured much more accurately than the ILC-I alone case, mainly 
because of the large VrVT-fusion cross section at ^/s = 1 TeV, see fig|3l The correlations 
between the error of c^^r^r and those of C2yyiY and 03^^^ are reduced because of the strong 
energy dependence of the contributions from the higher dimensional operators, as discussed 
in section IIVI In contrast, the reduction of the error of c-^^z is marginal and the correlation 
between c^^r^r and c^zz grows to 0.49, reflecting the dominance of the luminosity uncertainty 
A//2 

eq. fll78l) . mainly because of the strong energy dependence of their contributions; as discussed 
in section IVll The error of 03^^ is roughly half of the ILC-I result in eq. (11781) . partly because of 
the 2.3 times larger cross section of the no-tag {ee)H events, see fig|3]and eqs. (ll40p and (I14ip . 



0.005. The errors of c^zz 

factor of 3 smaller than the ILC-I results in 
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and partly because of the contributions from the double tag ee H events. 



B. Comparison with other papers 



First, we would like to compare our results with those of the reference [13|, as the present 
work is envisaged as an extension of this piece of work. The authors of ref. 13| analyzed 
the ZH production process for rriH = 120 GeV both with unpolarized and 90% polarized e~ 
beam in e"*" e~ collisions, using all the Z-boson decay modes. The study was based on the 
optimal observables method by allowing all the couplings to vary simultaneously. Although 
the formalism presented in section |V] of this report is significantly more compact than that 



of ref. [13|], we reproduce all the errors and their correlations for both unpolarized and 90% 
polarized e~ beams at ^/s = 250 GeV as given in eqs. (5.4) and (5.9), respectively, as well as 
eg. (5. 13) of ref.[3], for = 500 GeV with 90% polarized e~ beam. Note that the authors 



in ref. 13| considered the integrated luminosity 10 fb ^, and therefore the errors in eqs. (5.4), 



(5.9) and (5.13) should be multiplied by 1/ ylO for comparison with our results. It is worth 



noting here that the r lepton polarization and t 



le b jet charge identification with the efficiency 



40% and 20%, respectively, considered in ref. 13| has little impacts on the final results once the 



e beam polarization is available. The reduction of the beam polarization from 90% considered 
in ref. 13| to 80% in this report does lead to a slight increase in the error of the HZ'-f couplings 



by about 13%. 

Next we compare our results with tables II and III of ref.[l5|, where the authors estimated 
the 3(7 bounds on c^zz^ '^izzi '^rww ^^"^ '^2ww ( ^^z-, A6z/2, l^aw and A6vf/2, respectively, in 
their notation). This analysis was performed for unpolarized e"*" e~ collisions at ^/s = 500 GeV, 
with L = 500 fb~^ and rrifj = 120 GeV. It is worth noting that to derive these bounds they 
vary one coupling at a time and hence the inclusion of the 1% systematic error (accruing from 
luminosity uncertainty etc.) dominates their fiuctuation estimation for all the observed cross 
sections. Their results can be easily reproduced from ours simply by setting all the couplings to 
zero except for the one whose error is estimated, in the respective function for the specified 
process with an appropriate integrated luminosity, which should be corrected for the Higgs 
boson decay branching fraction. The systematic error of 1% is then added to the statistical 
error in quadrature. For instance, the 3 a limit \cizz\ — 0-034 [l^ is found from the double-tag 
eeH process excluding the Z —>■ e~ events. To compare, we set all the other couplings to zero 
in the function and find c^zz = ± 0.0092 for L = 500 fb"^ with B{H ^ hi) = 0.9. Now, 
adding the systematic error of 1% in quadrature we find c^zz = ^ 0.0105 and the 3 cr limit is 
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\(^\zz\ — 0.032. The limit \c2zz\ ^ 0.0022 15] is obtained from the observed ZH production 
cross section, where Z decays to muon pair and hght quarks (excluding b quark) pair have 
been considered. Under similar conditions we find the statistical error of Cg^^ = ± .00059. 
Since the luminosity error does not affect the measurement of Cg^^, this gives the 3 a limit of 
\c2zz\ < 0.0018, which shows an improvement by a factor of 1.3. We find that this difference 
is due to the use of the optimal observable, i.e., the differential distribution which is linear in 
C2zz'^ the quoted limit in [l5| is reproduced if we use only the C2zz effects on the total cross 
section. Likewise, we reproduce the bound jc^vFii^l' while we find a factor of two better bound 
on C2yy^r , whlch can be attributed to our use of optimal weight function. 



We also compare our results on c-^^^ and c^^^y couplings with those of the ref.[l7|. Since 
they studied possible constraint on the operator 0^2 {Oi in their notation), we present our 
comparison result in the next subsection; see the second footnote. 



C. Constraints on the dimension six operators 



The constraints on the effective HVV couplings given in eqs. fllTSp and (11791) should be 
expressed as those of 8 dimension-six operators of eq. ([2]) in order to compare the power of ILC 
precision measurements with that of the other experiments. All the effective HVV couplings 
are linear combinations of the coefficients /j/A^'s of these operators, as given in eq. f|T2|) . It 
is clear from eq.( fT2i) that the coefficients f^i, /^g, /<^4 cannot be determined uniquely from the 
HVV coupling measurements alone^^'. We therefore present constraints on two combinations 
of the three coefficients, f^^ and 3/^4 — 2f^2- Our results are hence for the 7 coefficients; /^^, 

/bw /w' /b' /wW' /bb ^^'^ 3/^4 — 2/^2- 

The combined ILC-1 results of eq. (ll78l) lead to the following constraints on the dimension-6 
operators, 

\ A ) ■I'pl 

V A J 
/ITeVNS r 

V A ) ^WW 

V A ) -IBB 



(^) (3/,4-2/, 



<f>2) 



= ± 


.091 


/ 1 










\ 


= ± 


.35 


-.40 


1 










= ± 


.051 


.37 


.23 


1 








= ± 


.084 


.23 


-.49 


.35 


1 






= ± 


.22 


.59 


-.91 


-.20 


.26 


1 




= ± 


.56 


-.49 


.98 


.22 


-.41 


-.96 


1 


= ± 


.36 


\-.052 


-.017 


.060 


.010 


.053 


-.032 1 / 



;i80) 



For instance, the measurement of the triple Higgs-boson coupling is necessary to constrain all the three 
operators. 
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We find that the three coefficients, f^^, and /g can be constrained rather accurately, with 
around 5% accuracy for A = 1 TeV, and rather independently of the other operators. The 
coefficient /^^ is measured accurately, because it contributes with the opposite sign to c^i^yy 
and Cizzy see eq. fll2al) and eq. (]12bp . In other words, f^^ is a measure of the difference between 
the HZZ and the HWW coupling strengths. On the other hand, the error of the 3/^2 — 2/^^ 
is rather large, since it measures the overall strengths of the HVV couplings. In fact, we 
find that it is only the error of 3/^2 — 2/^4 which is affected significantly by the luminosity 
uncertainty of A/ = 0.01. When we set A/ = 0, its error is reduced to 0.14™, but none of 
the errors of the other operators and their correlations are affected significantly. In addition, 
the three coefficients, /bw /ww /bb poorly constrained while their errors are strongly 
correlated, with the correlation matrix elements of —0.91, 0.98, —0.96, suggesting the presence 
of their linear combination, which can be measured accurately. We therefore present constraints 
on the 6 operators, after integrating out the contributions from 3/^2 — 2/^4 eq.( ll80p . in terms 
of the eigenvectors and their errors : 

- .096/^1 - .29/bw - -Se/w + -SI/b + .73/ww + •47/bb = ± -026 (181a) 

-.35/^1 - .18/bw + •84/w - -28^ + .21/ww + -IS/bb = ± -029 (181b) 

-.49/^, + .63/bw + -Orr/w + .47/b + .22/ww - •32/bb = ± -039 (181c) 

.59/^1 + .020/bw + •47/w + .65/b - .078/ww + -062/33 = ± .088 (181d) 

.52/^, + .49/bw + -OST/w - .43/b + .53/ww - •093/bb = ±.11 (181e) 

-.064/^, + .50/bw + •016/w - .052/b - .31/ww + -SI/bb = ± -69 (181f) 

As anticipated, we find 3 combinations of the 6 coefficients, whose errors are smaller than 
5% for A = 1 TeV. The worst constrained combination of eq. (1181fp has a much larger error 



of 68%, showing the poorly constrained combination of the three coefficients, /bw /ww 
/bb; which leads to their large errors and the strong correlations among themselves in eq. fll80p . 
It is worth noting here that the eigenvector of the most accurately measured combination in 
eq. fll81ap has a significant contribution from the constraint on C2^^, which is proportional to 



[% The authors of ref.[17| studied possible constraints on the dimension six operator 0^2 [Oi in their notation) 
which affects the HZZ and HWW couphngs by considering ZH production, WVF-fusion and double-tag 
eeH processes. They found the uncertainty of .005 for the coulpling 1c^^^ ~ 2c^-^-^ ( oi in their notation) 
at = 500 GeV for i = 1 ab"^ with B{H -> bb) ^ 0.9 and 80% fo-tagging efficiency. Our rcsuh of 
3/^4 — 2/^2 — ±.14 corresponds to 'ic^zz " '^'^iww ~ ±-0042 even though the total integrated luminosity 
of our analysis is 700 fb^^. This improvement is mainly due to the optimal observable method, but the 
luminosity uncertainty will limit our measurement. 
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/bw ~ /ww ~ /bb eq. fll2ep . In fact, if we drop the no-tag {ee)H events from the analysis, the 
eigenvector with smallest error becomes essentially that of eq.f ll81b| with a dominant term, 
which contribute to the c^zz^ C3Z7 ^zww couplings in eqs. p2g|)12ip . Also as expected, none 
of the results of eq. fllSip . neither the eigenvectors nor errors, are affected significantly by the 
luminosity uncertainty of A/ = 0.01. 

It is worth reporting here the importance of the e~ beam polarization to obtain the previous 
results. By setting |P| = 0, we find for the same lLC-1 integrated luminosities and their errors 
of A/ = 0.01, 



\ A ) 

v2 



( 



1 


TeV 




A 


1 


TeV 




A 


1 


TeV 




A 


1 


TeV 




A 


1 


TeV 




A 


1 


TeV 




A 


1 


TeV 


A 



BW 
W 



) /ww 
) /bb 



= ±0.17 




1 












= ±1.6 




.12 


1 










= ±0.22 




.45 


-.051 


1 








= ±1.2 




-.62 


-.089 


-.91 


1 






= ±0.61 




.41 


-.83 


.21 


-.24 


1 




= ±2.1 




-.031 


.986 


-.099 


.005 


-.91 


1 


= ±0.38 


v 


.13 


.30 


-.078 


-.009 


-.14 


.26 1 / 



(182) 



(^) (3/,4-2/,2) 

It is striking to find that all the errors except that of 3/^2 ~ 2/^4 are larger by more than 
a factor of 3 to 9 for the same luminosity. This is essentially because of the incapability to 
resolve the non-standard HZZ and HZ'-f couplings in the absence of beam polarization. On 
the other hand, we notice extremely strong correlations among the errors of /bw /ww /bb) 
and moderately strong correlation of -0.91 between the errors of fw and Jb- We therefore give 
the eigenvectors and their errors for the three most accurately measured combinations after 
3/^4 - 2/^2 is integrated out ; 



.24/^1- .43/bw± .27/w±.052/b±.66/ww± -50 /bb = ±-028 
.44/bw± -64 /w± .12 /b± .052 /ww - -30 /bb = ±-035 
,15/b- .12/ww+ -045 /bb = ±.083 



-•53/<^i 

.72/^,- .092 /bw± .66 /w± 



(183a) 
(183b) 
(183c) 



It is remarkable that the error of the most accurately measured combination in eq. fll83al) is 
not much different from that of eq. fll81al) . The reason is partly because both of them receive 
dominant contribution from the no-tag {ee)H events that measure the H'j'j coupling, which 
does not depend on the beam polarization. The second and the third combination of eq. fll83bp 
and eq. fll83cp have dominant contributions from the /^ and f^-^^. Except for the these three 
combinations, all the other eigenvectors have errors larger than 0.1 for A = 1 TeV. 
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Finally, our results for the combined ILC-I and ILC-II analysis, eg. (II 790 gives 
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.005 


-.005 


.029 


.002 


.028 


-.014 1 / 



;i84) 



i- 

When compared with the ILC-I only results of eq. fll80p . we find that the errors of and /g 
are reduced to 1/3, those of f^^, /ww ^^'^ /bb reduced to 2/3, while that of 3/^2 ~2/^4 
remains the same. The reduction of the errors in /g and is a result of the strong constraints 
on the c^iv^r, c^zz ^"^^ ^zz-y couplings at high energies. The strong correlations among the 
errors of /ww ^^"^ /bb remain unchanged, suggesting the persistent importance of the 
if 77 measurement via no-tag {ee)H events. The error of the combination 3/^2 ~ 2/^4 does not 
change, because it is dominated by the luminosity uncertainty. If we set A/ = 0, it reduces to 
±0.070. 

As in the ILC-I only case, we obtain the eigenvectors and their errors after integrating out 
the uncertainty in 3/^2 ~ 2/^4, and hence also over the luminosity error. We find. 



± 


.012 


(185a) 


± 


.015 


(185b) 


± 


.021 


(185c) 


± 


.038 


(185d) 


± 


.062 


(185e) 


±.43. 


(185f) 



-.12/^,± .001 /bw± .95 /w- .29 /b- .005 /^w" -028^6 
-.085/^1- .48/bw± -057 /w± .16/b± .66 /^w + -^^ Jbb 
-.27f^,+ .36/bw± .21 /w± -84 /b± .15/ww- -IT/bb 
.70/^,- .35/bw+ .22 /w± .42 /b - .39 /^w + -13 /bb 
.64 /^i± .53/bw± .043 /w- .14/b± .54 /^w" -063 /bb 
-.066 4i± .50/bw± .008 /w- .016 /b- .31 /^w ± -81 /bb 

We now find that 5 combinations out of 6 coefficients are constrained better than 5% for 
A = 1 TeV. We noticed that the most accurately measured combination in eq. (ll85ap is now 
dominated by refiecting the strong ^/s dependence of the c^ww^ '^szz ^^id c^^-y couplings; 
see eq.( 12g) to eq. fll2ip . The reduction of the error from that of the corresponding combination 
in eq. (]181bp is 60%. The second best constrained combination, eq. (]185bp . can be identified 
as the C2^^ combination, whose error is reduced by 40% from the ILC-I result of eq. (ll81ap . 
The coefficient of /b dominates the third accurately measured combination, eq. (ll85cp . The 
worst measured combination in eq. (]185fp is exactly the same as that of ILC-I only result in 



70 



eq. (1181f[) . while its error is reduced to about 2/3 reflecting a factor of 3 larger cross sections 
of the dominant WW- and ZZ-fusion processes at ^/s = 1 TeV as compared to those at 
= 500 GeV, see flgEl 



D. Comparison with the precision electroweak measurements 



A clear advantage of using the higher dimensional operators to parametrize possible new 
physics contribution is that we can compare the sensitivity and complementarity of any exper- 
iments, whether at high energies or low energies, in a model-independent manner. 

Although the HVV couplings can also be measured at the LHC in the Higgs-strahlung 
processes {WH and ZH production) and in the weak boson fusion processes, the expected 
sensitivity to the higher dimensional operators 33|, |3J] is not competitive with that expected 
at the ILC. On the other hand, the sensitivity of the precision measurements of the Z-boson 
and the H^-boson properties on the higher dimensional operators will remain competitive even 
in the ILC era. In this last subsection, we therefore compare our results with those of the 
present and future precision electroweak measurements. Although the results from LEP and 
SLC experiments have been flnalized 35|], both the mean values and the errors of the coefficients 
of the two operators, O^i and Obw in eq-©, will depend not only on nifj but also strongly on 
the continuously improving measurements of and m^r, and to a lesser extent on «s(?7i^)j;jg 
and a{'m\). We therefore present details of the dependences of the precision observables on 
these parameters. 

It is well known that the two operators O^i and Obw in eq.® contribute to the Z and 



ly-boson properties 0, l36|] via the oblique parameters 5* and T [23 



37 



38| 



(A5)np 

(Ar)NP 



"4 vr — /bw) 

1 „ 

J(j>i- 



186a) 
186b) 



2a A2 

Here (AS')np and (AT)np are the new physics contributions to the S and T parameters 



respectively. A 
ASz and ATz 



1 the Z boson parameters are parametrized in terms of the two parameters. 



24( 1 . which are related to the S and T parameters as 

A5 + Ai?z, 
AT + 1.49 Ai?z, 



^Sz 



(187a) 
;i87b) 



where ARz denotes the difference in the effective Z-boson coupling, QziQ"^) 22|, between = 
(where S and T parameters are defined) and = m\ (where the coupling is measured precisely 
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at LEP and SLC). By using the 13 data set of Z-pole parameters 35|, we obtain the following 
fit in terms of 4 parameters, AS'^, AT^, rrit and <ys{fnz)MS ! 

1. 



ASz = 0.037 + 0.0045 Xt - 0.037 Xs ± 0.105 

ATz = 0.043 + 0.0084 Xt - 0.065 ± 0.136 \ 0.90 1. 

/ X-t -+- X.. + z.^ \ ^ 

15.5 



A mm 



4.3 



0.17 



0.79 



188a 



;i88b) 



Here = (m^ — 172GeV)/ 3GeVandxs = (a^— 0.118)/ 0.003. The direct dependence of the 
fit comes from the Zhh vertex correction, and the dependence comes from T{Z ^ hadrons). 
In the above fit, AS*, AT and ARz are measured from their reference values in the SM at 



ra^ = 100 GeV, 



172 GeV, aJm 



zJms 



0.118 and Aa^^^ 



0.0277. Because the values 



of nifj, nit, «s and Aa^^^^ will be measeured precisely in the future, we parametrized their 
dependence as j2J, |39| , 



(A5)sM 
(AT)gM 



0.0963 Xh - 0.0224 xl + 0.0026 xl - 0.0014 Xt - 0.033 x^, 
-0.0432 Xh - 0.0539 xl + 0.0096 xl + 0.0367 Xt 
-0.0007 Xh Xt - 0.0033 



(189a) 

;i89b) 
(189c) 



{ARz)sM = 0.00838(1 -6-'"''). 

where, Xh = ln(m^/ 100 GeV) and Xa = (Aa^^^ - 0.0277)/ 0.0003. 

In addition to the Z-boson parameters, the VT-boson mass is also sensitive to the operators 
/^^ and /bw Their dependences can be parametrized as 

myy[GeV] = 80.318 - 0.288 (AS) + 0.418 (AT) + 0.337 (AU) - 0.0055 Xa, (190) 

where the SM contribution to the U parameter is 



(AU) 



SM 



0.2974 Xh - 0.0260 xl + 0.0772 xt + 0.0004 x 



t ) 



We note that none of the dimension-6 operators in eq.([2]) contribute to AU and ARz; 

(At/)NP = (Ai?z)NP = 0. 



(191) 



(192) 



latest estimates of m^^, m^, a 



By using the fit eg. (11881) of the LEP and SLC results [35| on the Z parameters, and the 

, Q and AaLd M 

m^;^,[GeV] = 80.403 ± 0.029, (193a) 
mjGeV] = 172.5 ± 2.3, (193b) 

(193c) 
(193d) 



C(s{mz)us = 0-118 ± 0.002, 



0.02768 ± 0.00022, 
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we find the constraints 



\ A ) 



BW 



0.048 ±0.028 
- 0.093 ±0.14 




(194) 



with Xmin/d-o.f. = 20.5/12, for = 120 GeV {xh = 0.1823). By comparing eq.ffTMD with 
eq.( ]180p for ILC-I and eq.( 1184p for the combined ILC-I and ILC-II analysis, we find that the low 
energy data constrain f^^ and better than the ILC. On the other hand, the two combination 
of the operators that are constrained by the low energy data 



.987 - .16 /bw 
.16 /^i± .987 /bw 



- 0.032 ± 0.015, 

- 0.10 ± 0.14. 



(195a) 
(195b) 



are quite orthogonal to the most precisely measured combinations at ILC-I, eq. (11801) . and those 
in the combined ILC-I and ILC-II analysis, eq. (11841) . For instance, if we combine eq. (ll84l) and 
eq. (11941) . we find 



( 1 TeV A^ n 

\ A ) J4>i 

V A ) JBW 

/1TW\2 n 
[ A J 

/1TW\2 r 

/'1TVV\2 r 

\ A J ^WW 

/1TW\2 r 

{ A I JBB 



- 0.048 ± .020 / 1 

- 0.093 ± .10 .68 1 

± .015 .27 .28 1 

± .025 -.12 -.17 .38 

± .065 -.34 -.65 -.40 

± .15 \ .59 .93 .36 



1 

-.25 
.009 



(196) 



-.86 1 / 



and now all the 6 operator coefficients are constrained rather independently, except for /bb 
which is still correlated with /b^ and /^w remarkable that even the errors of f^^ and 
/bw cire reduced significantly (~ 30%) by the ILC data. Because of the strong correlation 
among the errors of /bw /ww ^"^^ /bb ™- eg- (11^41) . the addition of the low-energy data (11941) 
leads to reduction of /ww ^^"^ /bb errors. 



We note here that the result given in eq. (ll96p does not take into account improvements 
in the measurements of m^r and which should certainly take place at the ILC-I, and also 
possible improvements in the measurements of asiji^z)uE ^'^Ld are expected. If we replace 
the present constraints of eq. (11931) by 



m,4,[GeV] 
mjGeV] 

as(?^z)MS 

^had 



80.403 ± 0.010, 
172.50 ±0.10, 
0.1180 ± 0.0010, 
0.02768 ± 0.00010, 



(197a) 
(197b) 
(197c) 
(197d) 
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without changing their mean values, the low energy constraints will become 

(1^)2^/^^ =_ 0.082 ±0.018 
(i^)'A/w =-0.21 ±0.12 




'198) 



Most importantly, the error of /^^ is reduced to about a half, because f^^ has rather strong 
dependence on rrit. 

In addition, if there are new measurements at GigaZ, we can measure the effective weak- 
mixing angle much more accurately. For instance, an estimate in ref. 4l|| gives 



sin^^^ = 0.23153 ±0.000013. 



(199) 



The effective mixing angle can also be parametrized as 2J, |39| 

sin^ 9^ = 0.23148 ± 0.00359(ASz) - 0.00241(Arz) ± 0.00011x„, (200) 



and the constraints on f,^ and will become 



0.097 ± 0.011 
0.32 ± 0.062 



1. 

0.95 1. 



(201) 



with xVd.o.f = 3O.9/I3M . Both errors of /^^^ and are reduce by a factor 2.5 and 2.3, 
respectively, from eq. fll94p . The eigenvectors and its errors are 



•98 f^, 
■17 /^i 



•17 /bw 
•98 /bw 



- 0.041 ± 0.0034, 

- 0.33 ± 0.062. 



(202a) 
(202b) 



Although the most accurately measured combination of the dimension-6 operator does not 
change from eg. 01951) . the error is reduced by a factor of 4.4. We note, however, that the reduc- 
tion of the errors are limited by the error of Aajj^^j assumed in eg. (11971 1). whose contribution 
to the uncertainty of the Xa term in eg. (12001) is 3 times larger than the error in eg. (ll99p . 
If we combine eg. ([MI) with ILC-I ± ILC-II result of eg.dHS]), we find 
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V A ; ^01 

/'lTeV\2 f 

\ A I -IbW 

( I TeV ^2 „ 
\~^^) 

\ A J JWW 

( \ TeV \2 f 

\ A ) ^BB 



0.097 ± 0.010 
0.32 ± 0.056 
± 0.014 
± 0.025 
± 0.053 
± 0.096 



/ 1 
.94 
.17 
-.088 
-.38 
V .75 



1 

.18 1 

-.093 .43 

-.42 -.36 

.80 .31 



1 

-.40 
.16 



M 1 



(203) 



(=*] Neither the low probabihty of the fit nor the \2>a evidence for (A/^;^,A/g^) in eq. (|20ip is our concern, since 
they are artifacts of our keeping the present mean values ea. (|197[) and eq. (|199p when reducing their errors. 
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Thanks to the precise measurements of f^^ and /bw errors of /gg and /ww reduced. 
The error of and /g are not affected much by the improved measurements of the weak- 
boson parameters, and the constraints of the type fll85ap ~ f ll85cl) will still give additional 
informations on new physics from the HVV coupling measurements. 

We hope that our report will be useful in studying the physics potential of the ILC project. 
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APPENDIX A: THE 3-BODY PHASE SPACE 



We parametrize the four momentum of the e^e — > f fH process, see eq. fll8p . in the labo- 
ratory frame as follows 

K = ^ (1,0,0, 1^, (Ala) 

K = ^ (1,0,0,-1^, (Alb) 

Pi = sin6'i cos 01, sin6'i sin0i, cos^i^, (Ale) 

P2 = fl, sin 6*2 cos 02, sin 6^2 sin 02, cos ^2)5 (Aid) 



2 

p'k = K + K- P'l - P'2, ' (Ale) 

= fl, sin^H, 0, cosO, (Alf) 



2 

where we ignore the masses of e^, / and /, and the Higgs-boson is produced in the xz-plane 

X\ sin 6*1 sin 0i + cos B2 cos 02 = 0. (A2) 
The Higgs-boson energy fraction and its velocity are 



2Ei 



xh = = 2 - xi - X2, (3h = xjl- ^m]j / sxj^. (A3) 
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We parametrize the / and / four momenta also in the rest frame of the // system, 
p*/ = !M (i^ sin r cos 0*, sin r sin 0*, cos^ 



rriff 
_ LL 

P2 — r, 



1, — sin 6** COS 0*, — sin^*sin</)*, — cos 6''' 



(A4) 



where m^j = (^1+^2)^, and (p* are the polar and azimuthal angle, respectively, with 
respect to the momentum direction of the // system in the laboratory frame. After making an 
appropriate boost to the laboratory frame with 

xhPh 1 



P 



7 



(A5) 



2-xh' ' v^r^2' 

and an appropriate rotation about the y-axis by 6^^, we find Pi and P2 in the laboratory frame 



m 



ff 



P2^ = ( - peg*), -Cg^Sg*C^* + Se^7(/5 - Cg*), -Sg*S^,, Cg^'j{p - Cg*) + Sg^Sg, 



(A6) 



Here we introduce a short-hand, sg = sin 6* and cg = cos 6*. Comparing eg. flAip with eg. flA6 
we find 



Xi 



ff 



l{l+PCg*), 



m 



X2 



ff 



7(1 - peg*), 



(A7a) 



„ Ce^l{P + Cg,)-Sg Sg*C^* Cg'j{p-Cg*)+S9Sg*C^* 

costal = — ^ ^ — , cos6'2 = — j-^ — ;r ; (A7b) 



1(1+ Peg*) 



sin 61 sin 1 



Sg*S^* 



sin 60 sin 1 



7(1 - /3cg* 
Sg*S^* 



l{l+Pcg*)' ^""'"^ 7(1 -/5C,.)' 

which satisfy eg. flA2p . We can now parametrize the 3-body phase space as 



d$3 = {2Tr)U^ {h + k2 - pi - P2 - Ph) 



d^Pi 



n (2vr)3 2E, 



d^PH 



(27r)3 2E^ 



(87r)2 
1287r3 



P 



m 



ff 



l\ d cos d m^^j d cos ^* ddf 



P 



dEjj d cos 9^ 



2-K All 
d cos 9* di 



1287r: 



■P 



dpTH dyn 



d cos 9* i 



s xi{l — xi — mjj / s) d cos 9i d(j)i d cos 92 d(j)2 

~ (JjUj ^ " 



1287r3 [l-xi(l-cos^i2)/2]2 



Att 



An 



(A7c) 



(A8a) 

(A8b) 

(A8c) 

(A8d) 
(A8e) 
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where 



P{a,b) = - 2(a + 6) + (a - 6)2, 

cos ^12 = 008 6*1 COS ^2 + sin ^1 sin 6*2 COS (</)2 — 



(A9a) 
(A9b) 



We use the phase-space parametrization (lA8bp in the analysis of the ZH production process, 
(IA8cP or flASdl) in the analysis of the v^VeH process where the undetectable Ve and Ve angles 
(cos6'*,0*) are integrated out, and flASeP in the analysis of the double-tag (eeH) events, following 
the prescription given in the Appendix A of ref. 27|. In the analysis of single-tag {e)eH events 
and no-tag {ee)H events, the electron mass should be kept in the integration of the forward 
scattering angles [27 1. 



APPENDIX B: MASSLESS FERMION CURRENTS 



In this appendix, we show the explicit form of the massless fermion currents which appear in 
the t-channel and s-channel gauge-boson exchange processes. The t-channel currents of eqs. fl47l) 
and (11131) for di = Xi = a and (72 = A2 = are 



J, 



e {cr,a) 



01 . ^1 



J 



'sxo sm — , COS — e 
' 2 2 



la COS — e 



-ia<p2 



2 



. ^2 

sm — 
2 



(Bla) 
(Bib) 



where cr/2 denotes the e~ helicities in eg. (IB lb ), and a/2 denotes the e"*" helicities in eq.f lBlb ). 
In the t-channel VT-exchange process, e"'"e^ Vf.Vf.H , only the a = — = — combination 
contributes, while in e^e~ — > e^e~H via t-channel Z and 7 exchange processes, both helicities 
contribute. The following 8 combinations of the contractions appear in the cross section with 
higher dimensional operators: 



3e-{--) 9tJ.u 3e+{+,+) 



Je-(- -) 9^lu Je+i--) 



3e-{+,+) 9tJ.u Je+(- ,-) 



^e-{+,+) 9^lu Je+(+,+) 



Je-(+,+)gi,.92/,Je+(- -) 



^(^1,^2), (B2a) 

^(^i,^2)m,^2,0), (B2b) 

= ^2)6^(^1, ^2,0), (B2c) 

= '^J'{9^,92)g'{9^,92A)■ (B2d) 
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Here, qi = ki—pi and q2 = k2—p2 are the transfer momenta in the t-channel, and 
The functions ^^(^1,^2), ^^(^1,^2,0), ^(^i,^2,0) and ^'(^i,^2,0) are 

9 6 

^{61 162) = 1s^x\X2 cos — sin — 

2 2 



7Y(6'i, 6^2, 



l_tan^cot^e^^ 
2 2 



2 - + cos^i) + xi sin cot yc"^"^ 
2 — X2(l — COS 612) + X2 sin ^2 tan -^^^'^ 

Q 

2 -xi(l + COS ^1) + xisinei cot ye*"^ 

Q 

2 — X2(l — COS 6'2) + X2 sin ^2 tan ^e*'''' 



>i -</>2- 

(B3a) 
(B3b) 

(B3c) 



(B3d) 



The currents that appear in the s-channel ZH production process are rather simple. The 
initial e^e~ annihilate currents are 



jI = V (k2, -|) 1^ P. u (ki, I) = v^(0, -a, -I, 0), 
where the e~ helicity is a /2. The final Z ^ ff decay currents are 



(B4) 



J2a' = U\Pl,^] Y Pa' V Ip2, 



y/s ^0, —a' COS 9* COS0* —i sin0*, —a' cos 9* sin0* +i cos0*, a' sin 6'*^ 



(B5a) 



in the // rest frame of eq. flA6l) . where the / helicity is a'/2. When contracted with the 
decaying Z-boson polarization vector 



e'^(A = ±) = -^(0,^1,-^,0) 
e/^(A = 0) = (0,0,0,1) 



(B6a) 
(B6b) 



we find 



e(±) ■ 32a' 

e{0)-j2a' 



2 



{l±a'ce*)e 



(B7a) 
(B7b) 
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The decay density matrix elements of eq. (l90|) are obtained from these equations. 
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